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We study the gravitational collapse of a magnetized neutron star using a novel numerical approach 
able to capture both the dynamics of the star and the behavior of the surrounding plasma. In this 
approach, a fully general relativistic magnetohydrodynamics implementation models the collapse 
of the star and provides appropriate boundary conditions to a force-free model which describes 
the stellar exterior. We validate this strategy by comparing with known results for the rotating 
monopole and aligned rotator solutions and then apply it to study both rotating and non-rotating 
stellar collapse scenarios, and contrast the behavior with what is obtained when employing the 
electrovacuum approximation outside the star. The non-rotating electrovacuum collapse is shown to 
agree qualitatively with a Newtonian model of the electromagnetic field outside a collapsing star. We 
illustrate and discuss a fundamental difference between the force-free and electrovacuum solutions, 
involving the appearance of large zones of electric-dominated field in the vacuum case. This provides 
a clear demonstration of how dissipative singularities appear generically in the non-linear time- 
evolution of force-free fiuids. In both the rotating and non-rotating cases, our simulations indicate 
that the collapse induces a strong electromagnetic transient, which leaves behind an uncharged, 
unmagnetized Kerr black hole. In the case of sub-millisecond rotation, the magnetic field experiences 
strong winding and the transient carries much more energy. This result has important implications 
for models of gamma-ray bursts. Even when the neutron star is surrounded by an accretion torus (as 
in binary merger and collapsar scenarios), a magnetosphere may emerge through a dynamo process 
operating in a surface shear layer. When this rapidly rotating magnet ar collapses to a black hole, the 
electromagnetic energy released can compete with the later output in a Blandford-Znajek jet. Much 
less electromagnetic energy is released by a massive magnetar that is (initially) gravitationally stable: 
its rotational energy is dissipated mainly by internal torques. A distinct plasmoid structure is seen 
in our non-rotating simulations, which will generate a radio transient with subluminal expansion, 
and greater synchrotron efficiency than is expected in shock models. Closely related phenomena 
appear to be at work in the giant flares of Galactic magnet ars. 



I. INTRODUCTION 

Understanding the gravitational collapse of a massive 
neutron star is of central importance for its connection 
to some of the most energetic astrophysical phenomena. 
Such an event may take place within a core-collapse su- 
pernova [1 , or in the late stage of a binary neutron star 
merger ^ ^3^ , and is widely believed to power some types 
of gamma-ray bursts. In recent years it has been re- 
alized that newly formed stellar mass black holes may 
be prodigious sources of electromagnetic energy, in addi- 
tion to driving strong kinetic outflows [U |4]-[6] . We are 
faced with the exciting possibility of probing the most ex- 
treme forms of gravitational collapse using coordinated 
measurements of electromagnetic transients and gravita- 
tional waves. Reflning models that make observational 
predictions for joint gravitational and electromagnetic 
radiation is critical in order to establish efficient obser- 
vation campaigns for both traditional astronomers and 
gravitational- wave astronomers. It is also a key step 
toward predicting the delay between peak electromag- 
netic and gravitational wave emission, the electromag- 
netic emission pattern (i.e. the beaming angle) and the 
electromagnetic and gravitational wave spectrum. 



Our focus here is on the evolution of the neutron star's 
magnetic fleld during its collapse to a black hole. We 
employ fully self-consistent relativistic calculations that 
follow the dense stellar material as well as the strong 
electromagnetic and gravitational flelds. A particularly 
intense electromagnetic transient is generated if the ini- 
tial magnetic field is very strong (~ 10^^ — 10^^ G), that 
is, if the star is a magnetar. Rapid rotation will enhance 
the energy of the transient, to a degree that can only be 
derived by a full time-evolution of the electromagnetic 
and gravitational fields. Rapid rotation also provides a 
context for generating strong magnetic fields: when the 
neutron star is accreting, the shear layer at its surface 
is a promising site for dynamo action. A distinct mag- 
netosphere will emerge and hold off the accretion fiow if 
even ~ 10~^ — 10~^ of the energy that is dissipated in the 
shear layer is converted to a poloidal magnetic field. Un- 
derstanding the evolution of the magnetic field around an 
isolated star is therefore potentially of key relevance for 
the collapsar and binary merger scenarios of gamma-ray 
bursts. 

Recent success in studying the behavior of plasmas 
around magnetized, spinning, stable neutron stars within 
fiat spacetime has been presented in Refs. [7]-[9]. To study 
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the collapse problem, one needs general relativity - to 
account for the role of spacetime curvature; and gen- 
eral relativistic magnet ohydro dynamics (GRMHD) - to 
determine the internal evolution of the star. The elec- 
tromagnetic (EM) phenomena outside the star can be 
approached in a variety of ways: through a full GRMHD 
calculation, which generally is very expensive given the 
low matter density; or more approximately using vacuum 
EM and force- free equations. In this paper, we compare 
all three approaches (restricting to ideal MHD in the case 
of GRMHD). 

Within general relativity, studies paying attention to 
neutron star collapse have been presented in the context 
of isolated stellar collapse (e.g. p!QHT4] ) and binary neu- 
tron star mergers (e.g. [15-20 ) with different degrees of 
realism. The role of electromagnetic fields are typically 
examined within general relativistic ideal MHD and stud- 
ies have been presented for single [21-26 and binary star 
systems [17, i27|. Such simulations have illustrated 
important details of the extreme dynamical behavior in- 
duced by the system which could trigger tremendously 
energetic phenomena. 

Newtonian and general relativistic simulations have 
shed light on how a binary neutron star merger can am- 
plify pulsar-strength magnetic fields by several orders of 
magnitude [18l I27H29] . The resulting hypermassive star 
is generically unstable to black hole formation, which 
opens up the tantalizing possibility that the increasing 
rotation rate and magnetic field strength would drive 
an intense electromagnetic outflow. To study such a 
scenario, our goal here is to combine the GRMHD ap- 
proach with a suitable description of the magnetically- 
dominated region outside a collapsing star, by applying 
the force-free approximation. Our approach is related 
to that of Ref. [30 , which matched a simplified, ana- 
lytic, relativistic solution for the interior of a collaps- 
ing star (dust-ball) with a numerical solution of the cou- 
pled Einstein-Maxwell equations in its exterior. Bound- 
ary conditions at the star's surface are provided by the 
analytic, ideal MHD solution. 

We go beyond this simplified scenario in three signifi- 
cant ways. First, our stars are evolved through collapse 
consistently by evolving the ideal GRMHD equations, 
and we can therefore consider in principle any kind of 
compact star. In particular, we have studied both rigidly 
rotating and non-rotating stars, whereas only the non- 
rotating case was considered in [30 . Second, the magne- 
tosphere is described within the force- free approximation, 
which, as argued in Ref. [ST, is a much more realistic 
model that, for instance, has been instrumental in un- 
derstanding pulsar spin down [8j. Third and last, our 
matching of the exterior solution with the interior one 
is dynamical. Thus the force-free solution can adapt to 
time-dependent fields sourced by the star. 

Our approach thus allows us to examine many interest- 
ing scenarios and we apply it here to study the behavior 
of collapsing, magnetized, compact stars (either rotating 
or not), the behavior of surrounding plasma (as described 



within the force-free approximation) and possible electro- 
magnetic radiation induced by the system. We compare 
our results with recent estimates in Ref. [32 for the non- 
rotating case which predict that the collapse process is 
smooth and that the magnetic field remains anchored 
to the star as a black hole forms, leaving a final black 
hole with a split-monopole field configuration. Our work, 
which follows the dynamics of the system, indicates that 
in both cases the stars radiate significant electromagnetic 
energy in which reconnect ion plays a crucial role, and 
that this radiation ceases shortly after the formation of 
a black hole. With force- free, such a black hole loses its 
electromagnetic hair within a dynamical time scale. 

We describe our hybrid approximation and the evo- 
lution equations in Sec. [Hj followed by a summary of 
the numerical techniques in Sec. HI Our choice of initial 
data is described in Sec. |IV[ and several tests of the hy- 
brid approach are presented in Sec. [V[ The new results 
for collapsing, magnetized stars are explained in detail 
in Sec. |Vl| while the astrophysical consequences are de- 
scribed in Sees. I VIII and I Villi We conclude in Sec. IIXI 
with some final comments. 



II. APPROACH 

In the presence of matter and electromagnetic fields, 
the Einstein equations must be suitably coupled to both 
the Maxwell and hydrodynamics equations. This cou- 
pling is achieved by considering the stress energy tensor 
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with contributions from matter and electromagnetic en- 
ergy given respectively by 
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A perfect fiuid with pressure P, energy density po, in- 
ternal energy e, and four-velocity describes the mat- 
ter state, and the Faraday tensor Fab describes the elec- 
tromagnetic field. The fiuid and electromagnetic com- 
ponents are directly coupled through Ohm's law, which 
closes the system of equations by defining the electric 
current 4- vector J a as a function of the other fields. A 
general relativistic expression can be obtained by consid- 
ering a multifiuid system of charged species [33 , leading 
to a fully non-linear propagation equation for the spatial 
component of the current. However, it usually suflices to 
consider a simplified version accounting for an algebraic 
relation between the current and the fields [34l 
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where cr''^ is the conductivity of the fiuid. The well- 
known scalar Ohm's law is recovered for a' 
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The equation of motion for the fluid and electromagnetic 
fleld are obtained from the conservation laws 

VaT-^ = 0; Va(po^^)=0; (5) 
VaF"^^ = ; Va*F'^^ = ; (6) 

which, together with the Einstein equations Gab = ^^Tah^ 
complete the system of equations governing the dynam- 
ics. 

Once the appropriate form of Ohm's law and the con- 
ductivity have been specifled, using for instance the al- 
gebraic relation of Eq. Q , the resulting equations typi- 
cally involve vastly different scales, rendering the imple- 
mentation of these equations quite costly from a com- 
putational point of view |114j . Fortunately, for speciflc 
regimes certain useful approximations capture the most 
relevant physics while bypassing the most strenuous dif- 
flculties. For our current purposes, involving magnetized 
neutron stars, the following relevant approximations can 
be deflned and employed in different regimes: 

• The ideal-MHD equations are obtained by requir- 
ing that the current remains flnite in the limit 
of inflnite conductivity, a ^ oo. This condi- 
tion also implies the vanishing of the electric fleld 
measured by an observer co-moving with the fluid 
{FabU^ = 0). This approximation is appropriate 
for the highly conducting matter expected in neu- 
tron stars. However, the numerical evolution of 
the ideal MHD equations typically fails in low den- 
sity regions where the inertia of the electromagnetic 
fleld is a few orders of magnitude larger than that 
of the fluid unless suflicient resolution is available. 
Such resolution requirements increase as the ratio 
of the electromagnetic to the fluid's inertia (or mag- 
netic to fluid's pressure) grows. Consequently, the 
approach becomes costlier in regions with decreas- 
ing physical relevance with respect to bulk mat- 
ter motion. Such a situation arises in particular 
in "vacuum" regions due to the standard practice 
of maintaining a density floor (a so called atmo- 
sphere) in regions of low density to exploit advanced 
numerical techniques for relativistic hydrodynam- 
ics. The atmosphere's density is much smaller than 
that inside the star, so this approach does not af- 
fect the relevant matter physics. However outside 
the star, the fluid inertia (or pressure) is typically 
much smaller than that of the electromagnetic fleld 
in magnetized cases and one generally encounters 
a large number of numerical difliculties. Different 
numerical strategies are often introduced to avoid 
them but such measures can limit one's ability to 
extract appropriate physics in these regions. 

• The force-free Maxwell equations are obtained by 
assuming the fluid's inertia is much smaller than 
that of the electromagnetic flelds. As a result 
FabJ^ = which in turn implies Fab^F"^^ = 
(^ E.B - 0) [3TJ[35]. This assumption therefore 



allows one to ignore the explicit time-evolution of 
the fluid as its dynamics are implicitly prescribed 
by the charge and current distribution. This ap- 
proximation is well suited to the dynamics of the 
low density, magnetically dominated plasma sur- 
rounding a compact object, but it cannot account 
for the physics in dense regions. 

• The vacuum Maxwell equations are trivially recov- 
ered assuming no coupling with matter (a = 0). 
This is a natural approximation in vacuum regions, 
far away from the compact objects of interest. 

As mentioned, computational costs are presently ma- 
jor obstacles to employ the general equations. It is thus 
highly desirable to deflne a scheme able to model all the 
relevant regimes in highly dynamical systems simultane- 
ously since none of these three approaches captures all 
the expected behavior. In particular we have in mind 
interacting binary neutron stars, black hole-neutron star 
binaries as well as the collapse of a magnetized star con- 
sidered here. All three systems are believed to play an 
important role in understanding gamma ray bursts and 
other energetic events driven by compact stellar-mass ob- 
jects. 

In what follows, we describe and apply a hybrid ap- 
proach which, while approximate, can account for the 
dynamical interaction of both gravitating and electro- 
magnetically driven fluids. Such an approach allows one 
to study the magnetosphere's behavior, in particular en- 
ergetics and fleld topology providing important clues for 
understanding relevant systems. To do so, we take into 
account that the electromagnetic inertia of a region with 
high conductivity (i.e. inside the star) would be orders of 
magnitude larger than that of the plasma region. Thus, 
we can ignore the back-reaction of the electromagnetic 
fleld in the plasma region onto that inside the star (i.e., 
a "passive" magnetosphere). We exploit this observation 
to deflne our approach, which employs both the ideal 
and force-free approximations suitably matched around 
the stellar surface. The matching procedure is such that 
the star's electromagnetic fleld provides the boundary 
conditions for the surrounding region treated with the 
force-free approximation, but the behavior of the mag- 
netospheric plasma on the high density stellar interior is 
ignored. 

Notice that this approach, in a sense, can be regarded 
as a natural extension of that adopted for studying 
pulsar-spin down by magnetosphere interactions [7, 8 to 
general relativistic dynamical contexts. In these works, 
only the behavior of the magnetospheric plasma is stud- 
ied and the star's influence is accounted for by boundary 
conditions. We here account for the possible stellar dy- 
namics and consequent influence on the boundary condi- 
tions deflned for the force-free approximation. 

In practical terms, we regard our system as described 
by two sets of electromagnetic flelds: (i) the ideal MHD 
flelds B^} and (ii) the force-free flelds {J^*, B*}, gov- 
erned by their respective equations. Both sets of flelds 
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are defined with respect to observers orthogonal to the 
spacelike hypersurfaces employed to foliate the space- 
time. Thus identifying the appropriate fields for match- 
ing is direct. The ideal MHD equations are evolved 
over the entire computational domain, as is customary. 
The matching region is determined by fluid density being 
some value above the vacuum region (though alternatives 
based on other physical quantities are obviously possi- 
ble). In particular, the force- free fields are evolved only 
in regions where po < p^^^^^. Dynamic boundary condi- 
tions are applied to the force-free fields on the surface at 
which po = p^^^^^ using the ideal MHD fields 



Kc = 



B' 



(7) 
(8) 



These conditions are applied in the spirit of penalty tech- 
niques [36 in which the equations are modified with driv- 
ing terms at boundary points in order to enforce the de- 
sired boundary condition. In particular, we extend the 
penalty technique across a number of grid points to ef- 
fect both the boundary conditions and to "turn-off" the 
evolution of the force-free fields within the higher density 
regions. Notice that this extension is not mathematically 
rigorous, but its usage will be justified by examining the 
solution in different test applications and comparing with 
the expected behavior. In practice we define a smooth 
kernel ^(^%^match) defined as 



where typically we adopt K ^ 0.001 /patmos and p"^^*^^ ^ 
200 - 2000 

Patmos? being Patmos the value for the density 
of the atmosphere. In the collapsing cases studied, this 
value is rescaled in time by the ratio of maximum density 
to the initial maximum density (see appendix). The val- 
ues found in Eqs. ( 7][8 ) along with the kernel Eq. ^ make 
their appearance in the equations determining {E\ B^} 
(together with "constraint cleaning" fields {^,^}, intro- 
duced to ensure constraints are well behaved through the 
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which are 
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The final term on the right-hand side of these equations 
is a penalty factor, which is introduced to impose inter- 
face conditions and ensure that a discrete energy norm 



is bounded. Details of this "penalty technique" are pre- 
sented in [36] and examples of applications in general 
relativity can be found in [38]. The above equations, to- 
gether with the Einstein and GRMHD equations are im- 
plemented as described in [25l [39^ „40^ to which we refer 
the reader for further details. 



III. NUMERICAL TECHNIQUES 

We adopt finite difference techniques on a regular 
Cartesian grid to solve the system. To ensure suffi- 
cient resolution in an efficient manner we employ adap- 
tive mesh refinement (AMR) via the HAD computational 
infrastructure that provides distributed, Berger-Oliger 
style AMR [H] \^ with full sub-cycling in time, to- 
gether with an improved treatment of artificial bound- 
aries [43 . The refinement regions are determined using 
truncation error estimation provided by a shadow hier- 
archy [44 which adapts dynamically to ensure the esti- 
mated error is bounded within a pre-specified tolerance. 
The spatial discretization of the geometry and force-free 
fields is performed using a fourth order accurate scheme 
satisfying the summation by parts rule, and High Res- 
olution Shock Capturing methods based on the HLLE 
flux formulae with PPM reconstruction are used to dis- 
cretize the fluid variables [151 US- The time-evolution 
is performed through the method of lines using a third 
order accurate Runge-Kutta integration scheme, which 
helps to ensure stability of the numerical implementa- 
tion [18]. We adopt a Courant parameter of A = 0.2 so 
that Ati = 0.2Ax^ on each refinement level /. On each 
level, one therefore ensures that the Courant-Friedrichs- 
Levy (CFL) condition dictated by the principal part of 
the equations is satisfied. 

To extract physical information, we monitor several 
quantities: (i) the matter variables and spacetime behav- 
ior, (ii) the electromagnetic field configuration and fluxes, 
and (iii) the electromagnetic Newman- Penrose (complex) 
radiative scalar (^2)- This scalar is computed by con- 
tracting the Maxwell tensor with a suitably defined null 
tetrad 



^2 = Fabn''m\ 



(14) 



and it accounts for the energy carried off by outgoing 
waves to infinity. The luminosity of the electromagnetic 
waves is 



-I 



dt 

lim / r^\^2\^dn . 



(15) 



Additionally we monitor the ratio of particular com- 
ponents of the Maxwell tensor 
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which, in the stationary, axisymmetric case, can be inter- 
preted as the rotation frequency of the electromagnetic 
field [35|. 



IV. INITIAL DATA 

Initial data for the hybrid equations involve the intrin- 
sic metric (gij) and extrinsic curvature (Kij) on a given 
hypersurface, as well as the magnetized fluid configura- 
tion in terms of its primitive variables (p, e, 5*). The 
initial data for the geometry and the fluid of rigidly rotat- 
ing neutron stars are provided by the LORENE package 
Magstar [46 , which adopts a polytropic equation of state 
P = Kp^ with r = 2 rescaled to = 100. Because the 
fluid inertia of a neutron star is many orders of magni- 
tude larger than its electromagnetic one, the magnetic 
field will have a negligible effect on both the geometry 
and the fluid structure, and so it can be specified freely. 
Unless noted otherwise, in our simulations we have cho- 
sen a dipolar structure for the initial magnetic field. The 
electric fields are set by assuming the ideal MHD con- 
dition Eq. ([7|), with zero fluid velocity in the exterior 
region. Additionally, we require initial data for the force- 
free fields {E^^B^}. Inside the star, they are defined to 
be exactly the same as their ideal MHD counterparts. 
Outside the star, the magnetic field is well defined by 
the dipolar solution, while the electric field is computed 
by assuming that the magnetosphere rotates rigidly with 
the star up to Re = 2Rs {Rs is the stellar radius), and 
imposing again the ideal MHD condition for the electric 
field. This configuration provides consistent data for the 
problem; however such data will not necessarily conform 
to the physical situation considered and so an unphysical 
early transient will be generated. 

V. TESTING THE APPROACH 

We first establish that the adopted procedure is in- 
deed able to capture correctly the dynamics of relevant 
systems. Since the ideal MHD approximation is self- 
consistently evolved throughout the computational do- 
main, our tests must address the behavior of the force- 
free fields. To do so, we first examine the convergence of 
our implementation and then illustrate that the approach 
provides the expected behavior by comparing with cer- 
tain recently studied cases. 

A. Non-rotating, magnetized star with a dipole 
magnetic field 

Adopting a stable and non-rotating stellar solution 
with mass M = l.GSM© and equatorial radius i?eq = 
8.62 km, from which we remove all initial pressure, we 
examine convergence of the force-free fields as the star 
collapses. This scheme contains "standard" sources of 
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FIG. 1: Non-rotating, collapsing solution. Results from FMR 
solutions at three different resolutions for the same nonrotat- 
ing collapsing star. Shown are the two differences between 
the resolutions. The difference in Bz along y = in the equa- 
torial plane at t = 0.06ms is first interpolated onto a uniform 
mesh for each solution. The difference between the medium 
and high resolutions is smaller than that between the low and 
medium resolutions, and this decrease as one increases the 
resolution indicates that the scheme is convergent. (Notice 
convergence in the central region is affected by initial data 
errors induced by depleting the pressure to induce the col- 
lapse.) 



error, such as: (i) truncation error of our finite differ- 
ence approximation of Maxwell's equations, (ii) errors 
associated with our application of the force-free condi- 
tions, (iii) the various numerical errors (truncation and 
constraint violations) associated with our GRMHD im- 
plementation as well as (iv) errors associated with our 
dynamic boundary condition that matches the force-free 
fields to the corresponding MHD fields inside the star. A 
detailed analysis of these errors is delicate and involved; 
however we illustrate that the fields obtained by this ap- 
proach converge to a unique solution with increasing grid 
resolution. 

We therefore evolve the collapse of a non-rotating 
star at three different resolutions with fixed mesh refine- 
ment (FMR). We subtract a high resolution run (with 
coarse level grid of 257^ points) and a medium resolution 
run (of dimension 193^) and do similarly for the medium 
and low resolutions (of dimension 129^). These two dif- 
ferences are plotted at t = 0.06ms (i.e., when the radius 
of the star starts to shrink) in Fig. [l] which shows that the 
difference decreases with increasing resolution, consistent 
with convergence. 



B. Magnetic monopole 

Next we consider a stationary force- free solution repre- 
senting a rotating neutron star with a monopole magnetic 
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FIG. 2: Monopole solution. Specific components of tfie mag- 
netic field displayed along the x = z = line at the initial 
and final times, together with the normalized density p/pc 
and the F function. Notice that the radial component has a 
smooth transition across the surface of the star, while there 
appears a toroidal component in the magnet osphere. 
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FIG. 3: Monopole solution. The rotation of the magnetic 
field lines (normalized with respect to its initial value) at dif- 
ferent times inside and outside the domain, separated by the 
kernel function F (continuous line). As time progresses, the 
rotation frequency approaches the constant value expected for 
the monopole solution. 



field [47], and compare to our results. To minimize effects 
due to oscillations of the star (induced by perturbations 
to the star induced by the discretization), the geometry 
and MHD fields are kept to their initial values and are not 
evolved. The initial data correspond to a rigidly rotating 
neutron star near the mass shedding limit, with a mass 
M = 1.84M0 and an equatorial radius i?eq = 12 km, 
rotating with a period of T = 0.886 ms. The magnetic 
field is given by = aBo{Rs/r)'^ , regularized conve- 
niently near the origin. As described earlier, the initial 
electric field satisfies the ideal MHD condition in a rigidly 
co-rotating magnetosphere which extends initially up to 
Re ~ 2 Rs. The evolution is performed in a cubic do- 
main of length L = 136 km ^ 14 i^^ with only two FMR 
grids and the star is placed at the origin. The maxi- 
mum resolution is Ax = 0.72 km, so that roughly 30 grid 
points cover the star. This resolution is sufficient for this 
test as both the geometry and matter variables are kept 
fixed and only the force- free equations are evolved until 
reaching a quasi-stationary configuration. 

The behaviour of different fields across the surface of 
the star after the solution has settled are displayed in 
Fig. [2j which plots radial profiles of: the kernel function 
F, the density po, and the non-trivial components of the 
magnetic field. As illustrated in the figure, the radial 
magnetic field only changes outside the star, where there 
also appears a toroidal component indicating the rotation 
of the magnetic fields in the magnetosphere. Overall the 
electromagnetic fields are seen to relax rather quickly to 
a state approaching the expected one. Another measure 
of the obtained solution is given by the (normalized) ro- 
tational frequency of the magnetic field in Fig. [3] As 
time progresses, this value approaches 1 as expected in a 
radially smooth way. 



C. Aligned rotator 

A particularly challenging test of our approach is the 
aligned rotator solution (see [8 ). The aligned rotator 
is a numerical solution of the force-free equations out- 
side a rotating surface representing a star that demon- 
strates closed field lines within the light cylinder (LC). 
This problem has been studied by a number of au- 
thors 111 m [48H5Q] working in flat space with a com- 
putational domain consisting only of the stellar exte- 
rior. In these works the star's influence is accounted for 
through suitable boundary conditions derived from the 
expected electromagnetic field at its surface. These ef- 
forts were specifically motivated to obtain the solution 
in the magnetosphere of isolated pulsars and, in particu- 
lar, to understand a possible spin-down mechanism that 
works even when the star's dipolar field is aligned with 
its angular momentum. 

As already mentioned, our approach differs in that 
we solve for the stellar interior, we work within curved 
space and our computational grid does not conform to 
the star's geometry [115^ . Because of these differences 
and also because our initial data for the force-free fields 
is only nonzero in the immediate neighborhood of the 
star, achieving the expected late-stage stationary solu- 
tion dynamically constitutes a demanding test. We note 
that the time and length scales required for such a test 
with respect to a realistic star are too computationally 
demanding if we are to capture the main physical as- 
pects of the solution (i.e. field topology, location of the 
light-cylinder, development of a current sheet, etc). In- 
stead of a realistic star, we adopt a rapidly rotating star 
so that the light cylinder is brought closer to the star - 
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FIG. 4: Aligned rotator. The fluid density and the mag- 
netic field lines on the x = plane at the four times t — 
(0, 1/3, 2/'^)T. Even before a complete rotational period, the 
solution exhibits the known properties of the aligned rotator 
solution. The light cylinder is located roughly at the expected 
position i?LC ~ 3.6 i?s (large tick marks indicate one stellar 
radius Rs). The intermediate plot illustrates the transient 
structure resulting as the LC forms and the initial data, which 
only extends to 2i?s, relaxes to fill the computational domain. 
These plots do not show the entire computational domain. 



where there is enough resolution to resolve it- and render 
the dynamical time scales shorter and easier to follow nu- 
merically. To achieve such rapid rotation, we resort to an 
unstable star but prevent the instability from disrupting 
the star by artificially freezing both fluid and geometry. 
Consequently we only evolve the force-free equations and 
compare the obtained configuration with the expected so- 
lution. 



We adopt a rotating star with a mass M = 1.7 Mq^ 
equatorial/polar radius of i?eq = 8.5/6.0 km and ro- 
tational period T = 0.64 ms. The light cylinder 
for this star is located at i?LC = c/^ = 31 km ^ 
3.6 i^s- The electromagnetic field is initially set simi- 
lar to that in the previous monopole test, but with a 
poloidal magnetic field given by a dipole outside the 
star. The evolution is performed in a cubic domain 
of L = 184 km ^ 22 Rs with four FMR levels with 
resolutions Ax = (0.25,0.5,1.0,2.0) km. The FMR 
hierarchy consists of centered cubes with side lengths 
L = (2.8,5.6, 11.2, 22)i?5, so that there are roughly 70 
points across the star in the equatorial plane. 

We evolve this star until the solution settles to a quasi- 
stationary solution. Fig. [4] illustrates the magnetic field 
topology at three representative times. In the last frame, 
we show the configuration to which the field settles, il- 
lustrating features predicted by known solutions for the 
aligned rotator. In particular, one can observe a region 
of closed field lines bounded at the expected radius, de- 
noting the light cylinder. 



VI. COLLAPSE OF A MAGNETIZED STAR 

Having tested our hybrid approximation to EM field 
evolution, we now focus on the collapse of a magnetized 
star. The star may or may not be rotating, the rotating 
case being of greater astrophysical relevance. In particu- 
lar, we are interested in studying how the magnetosphere 
responds, and the distribution of escaping electromag- 
netic radiation. 

Here it is essential to capture both the dynamics of 
the collapsing star and the surrounding magnetoplasma, 
since the two are tightly coupled. It is worth recalling 
that the introduction of rotation introduces a fundamen- 
tal difference in the radiative output of stationary (non- 
collapsing) stars as calculated in the electrovacuum and 
force- free limits: the spindown torque of an aligned rota- 
tor vanishes in electrovacuum but does not in the force- 
free case. Furthermore, this torque is strong and its mag- 
nitude varies only modestly with inclination angle [8 . 

Previous studies of the magnetosphere of a collaps- 
ing star are fairly limited. Some pioneering work in 
[5T] focused on the electrovacuum behavior of a non- 
rotating star. While this study included the effects of 
strong gravity, by neglecting plasma or rotation it was 
not able to capture the winding of the magnetic field dur- 
ing the collapse. The analytic study by [32 estimated the 
EM output during the collapse by applying the formula 
for the spindown luminosity of an equilibrium rotator, 
Lsd oc BIR^^VL^^ where Rs is the stellar radius, Bs the 
surface magnetic field, and Vt the spin frequency. How- 
ever, when the collapse time is shorter than the rota- 
tion period (as must be the case if the star is to avoid 
a rotational hang- up), the outer magnetosphere near the 
light cylinder is not able to follow the change in the sur- 
face magnetic field, and this equilibrium formula does 



not apply. Instead, twisting of the closed magnetic field 
lines is expected at a radius 

'max ctcoi/3 [116 , where 
tcoi = Rs/\Rs\ is the collapse time. The correspond- 
ing toroidal magnetic field is ~ (^^max/c)5(rmax), 
where 5(rmax) ^ {3Rs/ctco\)^Bs in a dipole geometry. 
The power injected into the magnetosphere is 



^sd 



3^ 

ct, 



col 



(17) 



which is larger than the equilibrium spindown power by 
a factor ~ {Qtco\/3)~'^ - 

To determine the dependence of L^d on Rg^ equation 
( p!7[ ) must be combined with the conservation of magnetic 
flux, BgR^ = constant; and the appropriate scaling be- 
tween Q and Rs {0. oc R~^ for self-similar collapse). The 
collapse does not, in general, follow a simple power-law 
relation between tcoi and i?s, but in the special case of 
pressureless collapse from a large radius {tcoi cx 
one obtains L^d oc R^^. It should be re-emphasized that 



equation (17) represents energy stored in the magneto- 
sphere, and so the maximum amount of energy that can 
only escape to infinity after the collapse is completed. 

During the last, relativistic stages of the collapse, ad- 
ditional physical effects arise. Spacetime curvature has 
a mixed effect on magnetic field winding in the magne- 
tosphere. On the one hand, the rotation frequency is re- 
duced with respect to the Newtonian value; on the other 
hand, torsional Alfven waves in the force-free magneto- 
sphere slow down significantly as the horizon approaches 
the surface of the star. 

The net result, as we show, is that the magnetic field 
becomes strongly wound up if the star is rotating close 
to breakup before the collapse, and develops a/M > 0.5 
after the collapse. To clarify the influence of plasma and 
rotation, and to make useful comparisons, we have made 
a parallel set of runs in the electrovacuum approximation. 
These are obtained straightforwardly by setting J* = 
and not enforcing the conditions {E.B = 0, < \B\} 
in Eqs. ([l0][l3|. 
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A. Non-rotating stellar collapse 
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FIG. 5: Non-rotating, unstable star (force- free) . The fluid 
density (colors in the central region), the magnetic field lines 
(blue) and the EM radiation flux density (red) at times 
t = (-1, -0.2, 0.02, 0.05, 0.1, 0.2) ms (from left to right, top 
to bottom). Here t — denotes the time that the horizon 
appears. The magnetic field lines are dragged by the star 
during the collapse, producing Alfven waves in the magne- 
tosphere that carry a small fraction of the magnetospheric 
energy and stretch the magnetic field lines near the equato- 
rial plane. Most of the EM energy falls into the black hole 
(the solid white, central rectangle denotes the excision region 
for the singularity). 



The collapse of a magnetized, but non-rotating, neu- 
tron star was studied in [30], where a magnetic field 
frozen into the stellar surface was matched to an exte- 
rior, vacuum solution of Maxwell's equations. This cal- 
culation followed the transition to a black hole, and the 
ringdown of the EM field threading the horizon. 

In reality, the medium outside the star is an excellent 
electrical conductor. As is argued in [3T1|35], a combi- 
nation of large voltages and strong gravitational fields 
will trigger runaway pair creation near compact objects 
[3T| l35] . In a strongly dynamic situation, the number 
of particles generated can be enhanced even further, e.g. 
by a Kolmogorov-like transfer of energy from large-scale 
waves to internal plasma heat [52 . The remnant of a bi- 
nary merger or a rapidly rotating stellar core collapse is 



also a strong neutrino source and, therefore, its magneto- 
sphere will be filled with a much denser baryonic plasma 
compared with pulsar magnetospheres |53l |54]. But in 
most cases, the plasma around compact objects is mag- 
netically dominated and its dynamics is nearly force free. 

Therefore, we expect significant differences in the time 
evolution between the electrovacuum exterior assumed 
by [30 and a force- free magnetosphere. To anchor the 
external force-free (or electrovacuum) solution, we adopt 
a marginally unstable, non-rotating star with a mass 
Ms = 1.63M0 and a radius Rs = 8.62 km. The nu- 
merical domain extends up to L = 16 Rs and contains 
three centered FMR grids of sequentially half sizes (and 
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hence twice better resolved). The highest resolution grid 
has Ax 0.18 km. 

The star collapses to a black hole in ~ 1 ms [25] . We 
set t = at the onset of an apparent horizon and dis- 
play the magnetic field and the density at various times 
in Fig. |5j This force-free solution has several salient fea- 
tures. There is an early transient in which the dipole 
magnetic field relaxes to a solution consistent with the 
physical configuration. Subsequently, as the star col- 
lapses, the magnetic field lines are gradually stretched 
along the equatorial plane. After 1.0 ms an apparent 
horizon appears in the interior of the star, which grows 
as it swallows all the remaining fiuid in ~ 0.15 ms. As the 
outer layers of the star are accreted by the black hole, the 
stretched magnetic field lines near the equatorial plane 
reconnect and form closed loops that carry away electro- 
magnetic energy and magnetic fiux (note the field loops 
in the last frame of Fig. [5| . 

The EM field evolution shows qualitative differences 
in the force- free and electrovacuum runs, especially af- 
ter the black hole forms (Fig. [6|. Starting from a com- 
mon dipole structure at the beginning of the collapse, 
the magnetic field in the force-free case becomes radi- 
ally stretched near the magnetic equator, and maintains 
a consistent sign. In the electrovacuum evolution, the 
stellar magnetic field disconnects more readily from the 
exterior. Changes in the connectivity of the magnetic 
field follow the appearance of zones where E'^ > B'^ ^ as 
is demonstrated in Figs. [lO] and [12] 

For a more quantitative discussion, we consider in 
Fig. [T] the electric and magnetic fiuxes, computed over 
a surface located at r = 1.5i?s, as a function of time. 
As expected, both the total (signed) fluxes remain small 
throughout the simulation, indicating that essentially no 
spurious magnetic/electric charges are created during the 
collapse. The unsigned magnetic flux decays exponen- 
tially in both the electrovacuum and force-free cases. In- 
terestingly, the decay rate in the former case matches 
the / = l,m = quasi- normal mode for electromagnetic 
perturbations of a Schwarzschild black hole with mass 
M = 1.63M0 [55 . 

The decay rate of the magnetic flux in the force-free 
run is roughly twice the electrovacuum result. One could 
translate the measured e-folding time Ie into an "effec- 
tive reconnection speed," Vrec ^ GMc~^tE~^ by observ- 
ing that reconnection takes place within r o:^ (1 2)rH 
(with rn the horizon radius). Our results indicate Kec — 
0.14c. 

This result may appear paradoxical at flrst sight, and 
deserves some comment. Electrovacuum magnetic flelds 
effectively reconnect (by converting to electric-dominated 
flelds) at the speed of light, as we discuss in more de- 
tail in Sec. |VI A 4[ Why then is the decay faster in the 
force-free case? The answer appears to reside in the self- 
inductance of the black hole. As magnetic fleld lines 
reconnect through the equator, they generate a strong 
toroidal EMF. In a vacuum, this EMF sources a mag- 
netic flux of the opposing sign. When conducting matter 
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FIG. 6: Non-rotating, unstable star (electrovacuum). The 
fluid density (colors in the central region) and the magnetic 
field lines (blue) at t = (0.03, 0.06, 0.12, 0.125, 0.13, 0.15) ms. 
The magnetic field lines are dragged by the star during the 
collapse, stretching the magnetic field lines near the equato- 
rial plane until that the fluid is swallowed by the black hole. 
Afterward, the EM dynamics is mostly described by the Quasi 
Normal Modes of the system. 



is present, there is no oscillation in the sign of the flux 
threading the hole. Rather, reconnection is a monotonic 
process and, after the formation of an x-point, the mag- 
netic fleld lines interior to the x-point fall through the 
horizon. 



1. Electromagnetic Output 

In order to estimate the efliciency with which EM en- 
ergy is radiated to inflnity, we have scaled the time- 
integrated luminosity to the peak electromagnetic energy 
contained in the magnet osphere. This peak energy is 
reached approximately at the formation of the apparent 
horizon. In the absence of rotation, it is Cpeak ^ 2 times 
the initial dipole energy, as determined by the conser- 
vation of magnetic flux. The measured radiated energy 
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is, in turn, a fraction e^ad of Cpeak-Edipoie.o- Numerically, 
i?dipoie,o « (27r/3)(i?2„iJ87r)i?3 = {l/12)Bl^,^Rl and 
one finds £^dipoie,o = 1-4 x B'^^^^^^^eig for a (initial) 
polar field 5poie and radius Rg ^ 12 km. Hence 



£^rad ~ 1.4 X 10"^^ Cpeak ^rad ^pole,15 ^rg. 



(18) 



The radiative efficiency is very small in the force-free 
case, erad = 0.008; the rest of the peak energy is swal- 
lowed by the black hole. This is illustrated in Fig. |8) 
which shows how erad grows with time in both the force- 
free and electrovacuum runs. The radiated energy is 

from Eq. (18), expressed in 



^rad « lO^^B^^j^ ^gergs 
terms of magnetar-strength magnetic fields. Most of this 
is radiated in a short interval « 1 ms surrounding the 
collapse, with an average luminosity L ^ 10^^ ^p^^^ ^5 
erg s~-^. 

The radiative efficiency measured at r = l.bRs is an or- 
der of magnitude higher in the electrovacuum case. This 
indicates that a larger proportion of the electromagnetic 
energy falls into the black hole in the force-free run, in- 
stead of escaping to infinity. An important check of this 
result is to measure the dissipation at current sheets in 
both the force-free and electrovacuum simulations. In 
both cases, we find that the integral of E.J is a small 
fraction of the total EM energy in the magnetosphere. 
(The constraint E'^ < is maintained in the force-free 
case by applying a small enhanced resistivity; hence the 
reduction in electric field energy appears as E.J dissipa- 
tion.) 

An important feature of this radiation in both cases is 
its predominantly dipolar structure in energy flux (i.e., 
L oc sin^(^)). In the force- free case, energy is radiated in 
a rather continuous manner and it propagates outwards 
with a velocity v = 0.89 c, equivalent to a mild Lorentz 
factor of = 2.2. For the electrovacuum case the energy 
is radiated mainly in two long bursts instead of the sev- 
eral, and with shorter periods, resulting in the force-free 
case. 



2. Late Force-Free Evolution and Magnetic Reconnection 

At later times, when the fluid has completely fallen 
into the black hole, the field lines that were dragged to- 
ward the horizon reconnect near the equatorial plane in 
a few sequential bursts, expelling most of the remaining 
magnetic flux. This behavior appears to be associated 
with the formation of x-type singularities in the magnetic 
field. Since the field is stretched radially and then re- 
connects near the horizon, the resulting electromagnetic 
pulse bunches up in the radial direction, which explains 
the structure seen in the last panel of Fig. [5] 

X-point reconnection appears to happen easily in force- 
free plasmas; in ohmic plasmas it is associated with in- 
homogeneities in the electrical resistivity (e.g see [56 ). 
In the absence of a detailed microphysical model for the 
resistivity, fine details such as these should be treated 
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FIG. 7: Non-rotating, unstable star. The absolute value of 
the magnetic flux in the electrovacuum case as a function 
of time, and the electric and magnetic fluxes in the force-free 
case. These quantities are integrated over central spheres with 
radius r = 1.5i?s and normalized with respect to the initial 
integral of |$s(t = 0)|. The total signed fluxes remain very 
small throughout the simulation, and the unsigned magnetic 
flux decreases as the black hole swallows all the matter which 
anchors the magnetic field. 
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FIG. 8: Non-rotating, unstable star. Time integral of the elec- 
tromagnetic luminosity, normalized with respect to the peak 
EM energy of the magnetosphere (i.e., around the formation 
of the black hole), in both the force-free and electrovacuum 
simulations. In the electrovacuum case, the net EM output 
is ^ 8%, ten times larger than that radiated in the force- free 
case. 



with caution, and the calculation should be viewed as 
illustrative. 



3. Vacuum versus Force-Free EM Field Evolution in 
Axisymmetric, Non-Rotating Collapse 

Some important features of the vacuum evolution of 
the electromagnetic field around a collapsing star can be 
understood by neglecting the effects of spacetime curva- 
ture, and by considering a simplified trajectory for the 
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surface of the star. If the cohapse starts at a finite time, 
then an initially potential magnetic field evolves into a 
hybrid structure that consists of an inner potential mag- 
netic field that matches the surface boundary condition 
as determined by the conservation of magnetic flux; and a 
transient electromagnetic wave that propagates into the 
original field structure. 

Of especial interest is the appearance of zones within 
this wave structure that are dominated by the electric 
field. If the collapse continues for a long time (the fi- 
nal stellar radius Rs is small compared with the initial 
radius), then this zone where > extends over a 
wide range of radius. Since a realistic magnetosphere 
may contain enough free charges to limit the growth of 
E • this provides a nice example of how the nearly 
force-free evolution of an electromagnetic field can lead 
to strong dissipation. 

In the absence of rotation, a spherically symmetric 
collapse implies a radial fiuid velocity inside the star. 
As long as the stellar surface contracts with a speed 
—Rs <^ c, the magnetic field near it is approximately 
potential. For a pure multipole of order ^, 



B{r) = Bs{t)Rs{t)W 



(19) 



(with Bs the magnetic field at the star's surface). A 
toroidal electric field E^j^ = —{Rs/c)Bo is present at the 
surface of the star, assuming its interior to be perfectly 
conducting. The junction condition at the surface en- 
sures the continuity of E^jy. The surface magnetic field 
therefore increases in accordance with simple fiux con- 
servation. 



dBs 

dt 



dBs 
dr 



Rs 



Rs 



(20) 



At a fixed radius r, the magnetic field grows weaker: in 
the case of a simple dipole, the stellar magnetic moment 
scales as /j.{Rs) = Bs{Rs)Rl = fio{Rs/ Rso)- 

If the star were to reach infinite density at a finite time 
tcoi, then a strong toroidal electric field would develop at 
r > c(tcoi — t). The inner potential zone would shrink 
along with the star as t approaches tcoi- 

The external electromagnetic field can then be ob- 
tained by rescaling the radius, r ^ = r/Rs{t)^ and 
transforming derivatives according to dtXir^t) [dt + 
{Rs/Rs)d^]X{£,^t). It is simplest to solve for the vector 
potential A^, from which the poloidal magnetic field and 
toroidal electric field are derived. The boundary condi- 
tion at the surface of the star is 



Rs{t)A4Rs{t),0] 



Bs{t)Rl{t) dPj ^ RsoA^ojO) dPj 
£ dO £ 



£ dO 

(21) 

From Eq. (20), RsoA(f)Q{6) is constant. Substituting 
A(f){r^6) = RsoA(f)o{6)g{(^,t) into the wave equation 



r 



(22) 



and adopting a collapse law Rs{t) cx (tcoi — t)^ (here 
a = 2/3 for pressureless collapse from a large radius), we 
find 



dig + -drg 

a 



c 



dig 



£{£ + !) 



(23) 



Here t = J adt/{tcoi — t) is a dimensionless time coordi- 
nate; hence Rs{t) = Rsoe~^ . The electromagnetic field 
is constructed from the solution to Eq. ( 23 ) using 



Be{r,t) 



^ — de (smM^o) ; 

sm^ 

RsoA(f)0 

-^iid^9)^=r/R/. 



Rs RsoA(i)0 



i=r/R. 



■ (24) 



One recovers the usual potential solution ^^(r, ^) = 
A4Rs{t),0]{r/Rs)-^^^^^ where ^ < c/\Rs\. A self- 
similar solution is also available in the case of collapse 
at a uniform speed, a = 1, if the collapse starts at a very 
large initial radius. Then one can take drg = and the 
electromagnetic field is a function only of r/Rs. In this 
case, the magnetic field can retain a dipolar form out to 
large distances r c{tco\ — t) from the star. Restricting 
to ^ = 1 gives 



Beirut) 



Rs 2/iocos6> 

Rs /io sin 
r 



c(tcoi - 1) 



Be{r,t). 



(25) 



Although the magnetic field is identical to that 
sourced by a stationary dipole of magnitude fi{Rs) = 
jiio{Rs/ Rso)^ the electric field energy dominates outside 
a distance ~ {c/\Rs\)Rs{t) from the star. The sim- 
ilarity solution is accurate out to the larger distance 
~ {c/\Rs\) RsO^ where Rso is the stellar radius at the 
beginning of the collapse. 

To understand how this inner solution for the electro- 
magnetic field matches onto the initial potential mag- 
netic field that was present prior to the collapse, or to 
consider cases other than constant Rs^ one must calcu- 
late the full time-dependent solution to ( [23| ). We have 
done this by evolving g{^^r), drg{£,^r) and 9^^(^,r) us- 
ing a centered discretization of Eq. (23) and employing 
a small Kreiss-Oliger dissipation value O(10~^) to damp 
high-frequency noise in each of these variables. The re- 
sulting electromagnetic field profile is plotted in Fig. [9) 
with distance normalized to the initial radius of the star. 
One observes in Fig. 10 the emergence of an extended 



zone with E"^ > B^, as is expected from the similar- 



ity solution of Eq. (25). The magnetic field is relatively 
stronger in the polar regions, since (x sinO ^ at 
small 0. 
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FIG. 9: Evolving vacuum EM field around a star collaps- 
ing to half its initial size (uniform dRs/dt). Top panel: the 
inner magnetic field tracks the instantaneous dipole of the col- 
lapsing star; at a fixed radius, B oc ii{Rs) oc Rs/Rso- Time 
progresses top to bottom. Bottom panel: The zone of rising 
E(f) closely follows the similarity solution of Eq. (25). Time 
progresses left to right. (Two colours are employed for clar- 
ity.) 



Eventually the amplitude of the outgoing wave distur- 
bance becomes large enough that Bg changes sign. After 
this happens, the magnetic field remains dipolar inside 
the radius ~ {c/\Rs\)Rso^ but disconnects from the ex- 
ternal zone of undisturbed potential field. This transition 



is illustrated in Fig [TT] using two snapshots corresponding 
to the two most extended field profiles in Fig. [9j 

The first part of the collapse leads to a re- arrangement 
of the magnetic field outside the star, while electromag- 
netic energy flows inward: the Poynting flux. 



Sr = —E^BqC '. 



(26) 



is negative inside r ~ {c/\Rs\)Rso- After the horizon 
forms and reaches the surface of the star, the compres- 
sion of the magnetic field stops, and some of the trapped 



FIG. 10: Relative strength of vacuum electric and magnetic 
fields around collapsing star (uniform dRs/dt). Different lines 
illustrate the obtained behavior for different values of (from 
10"", top, to 90"^, bottom). The zone where E'^ > B'^ is more 
extended near the magnetic equator. 



magnetic field can be radiated to infinity. The angular 
distribution of this radiated energy reflects the symme- 
try of the initial field. In the case of a dipole, the energy 
is carried away by closed magnetic loops and is mainly 
channeled through the magnetic equator. 

We note that the appearance of zones with E'^ > B^ in 
the electrovacuum solution implies a fundamental differ- 
ence with the alternative force-free solution. In this case, 
E-B = throughout the collapse, so consistent force- free 
evolution is obtained only by removing energy from the 
electric field. From an MHD perspective, such a tran- 
sition signals the appearance of nearly luminal plasma 
motions, where the inertia of even a small residue of en- 
trained matter can become important. 

This simple example shows that there can be profound 
differences in the macroscopic structure of the electro- 
magnetic fields when even a small amount of conduct- 
ing matter is present. In the complete absence of free 
charges, macroscopic zones where E^ > B^ are the natu- 
ral consequence of the time-evolution of an EM field that, 
initially, is purely magnetic. 



4. E'^ > B'^ in the Vector QNM of a Black Hole. 

The quasi- normal mode (QNM) behavior observed 
shows two remarkable features: an oscillation in the sign 
of the magnetic flux threading each hemisphere, and the 
appearance of an equatorial zone with E^ > B^ . These 
two features can be related to each other using a simple 
planar analogy. 

Consider an initial magnetic field configuration B = 
Bqx {—Bqx) for y > (< 0). This is familiar from stud- 
ies of conducting fluids, where in the case of uniform 
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FIG. 11: Connectivity of the magnetic field around a col- 
lapsing star (uniform dRs/dt). The dipolar magnetic field 
that is anchored in the star disconnects from the surrounding, 
undisturbed, potential field when the amplitude of the output 
wave becomes large enough that Be reverses sign. Note that 
> (the red sections of the field lines) in the outer part 
of the inner zone with a dipolar magnetic field line profile. 
These two frames display the solution at the two latest times 
shown in Fig. [9] in which Be becomes negative. 



resistivity one finds a slow flow of matter to a thin cur- 
rent sheet at = 0. The sheet thickness is limited by 
the flow of matter along the sheet to large \x\ (Sweet- 
Parker reconnection). But now we are interested in the 
case where the conducting matter is absent, and the elec- 
tromagnetic field evolves according to the vacuum wave 
equation. One finds, instead, a growing zone of pure 
electric field E = —Bqz for \y\ < ct, which maintains a 
uniform sign across the initial magnetic null surface. A 
pure magnetic field is converted to a pure electric field. 

In the case of EM fields localized around a black hole, 
this conversion of magnetic to electric fields occurs at the 
magnetic equator, where the poloidal field lines merge to- 
gether. The toroidal electric field that is created sources 
a poloidal magnetic fiux through the horizon of the op- 
posing sign to the pre-existing flux. In this way, a con- 
tinuing interconversion of magnetic and electric fields can 
be maintained. 
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FIG. 12: Non-rotating, unstable star (electro vacuum). 
Electrically-dominated regions {E^ > B^) are marked in 
color, at times t = (0.04, 0.13, 0.22, 0.31) ms. Here cur- 
rents arise in the force-free case, and < B^ is maintained 
through gradual dissipation, resulting in a fundamental dif- 
ferent with the electrovacuum evolution. As described by the 
simple Newtonian model, these regions form near the equato- 
rial plane and close to the collapsing star. As time progresses, 
they propagate outward in bursts. The grey zone in the center 
represents the star. 



B. Rotating stellar collapse 

The collapse of a rotating, magnetized star produces an 
interesting generalization of the relativistic wind problem 
for a stationary star (e.g [8l[31]). As in the case of spheri- 
cal collapse, qualitatively new effects are introduced after 
the formation of a horizon. We have chosen an unstable, 
rotating model star with a mass M = 1.84M0 and equa- 
torial radius Rs = 10.6 km. The star rotates with a 
period T = 0.78ms, so that the light cylinder is initially 
located at i?LC = 37 km ^ 3.5 Rs. The numerical do- 
main and resolution are identical to those employed in 
the non-rotating case. 

To remove unphysical transients we evolve the force- 
free equations for a couple of periods with both geom- 
etry and matter fixed, and afterwards all equations are 
evolved. This approach ensures the force-free fields relax 
to a configuration consistent with the physical scenario 
considered. The expected field configuration emerges 
during this startup phase, with a closed, corotating mag- 
netosphere extending out to the light cylinder (LC). The 
Goldreich- Julian (G-J) current structure is present, with 
an outflow along the polar field lines balanced by a return 
flow through an equatorial current sheet. 

The different stages of the collapse are represented in 
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FIG. 13: Rotating, unstable star (force-free). Magnetic field 
configuration (blue lines) and fluid density (marked in red) 
at times t = (-0.47, -0.17, -0.01, 0.12) ms. As the collapse 
proceeds the increasing spin rate of the star pulls the external 
magnetic field in the toroidal direction. 



4- 




FIG. 14: Rotating, unstable star (force- free) . Star's angu- 
lar rotational velocity, measured at the equator, during the 
collapse. 



Fig. [13) while that the angular velocity of the star is dis- 
played in Fig.[l4) As the star contracts, and its rotational 
frequency increases, the instantaneous LC approaches the 
star [1 17 . Differential rotation develops in the magneto- 
sphere, due to the lack of causal contact between the 
star and the LC, and the magnetic field is wound in 
the toroidal direction. Furthermore, the deepening of 
the gravitational potential forces significant changes in 
the magnetic field profile, by pulling the field lines more 
tightly toward the star. 

The poloidal magnetic field strengthens due to flux 
freezing in the star, just as in the non-rotating case, but 
now most of the EM energy is in the toroidal component. 
Near the poles, the field lines twist around, generating 
a cone-like structure. In general, the magnetic field pre- 
serves a stretched dipolar topology for a longer time than 
in the non-rotating case, up to the point that all the fluid 
is swallowed by the black hole. 

As the black hole forms inside the star, and the fluid 
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FIG. 15: Rotating, unstable star (force-free). Charge 
density (blue: positive and red: negative) is plotted at 
t = (-0.17, -0.01, 0.12, 0.35) ms. The central fluid region 
marked in black disappears as the horizon emerges. 



falls inward, many of the features in the EM field evo- 
lution are preserved from the non-rotating runs. Snap- 
shots of the field profile, charge density and Poynting 
flux, taken at various times close to horizon formation, 
are displayed in Figs. [15] and [16) Just as in the non- 
rotating case, the topology of the magnetic field lines 
changes dramatically: the y-point structure in the mag- 
netic field disappears from the equatorial regions, and a 
current sheet extends inward to the horizon. This current 
sheet is subject to spasmodic episodes of reconnection, 
the details of which may depend on the prescription for 
the electric resistivity and its variation with radius. 

In less than a millisecond, the magnetic flux thread- 
ing the horizon has almost completely vanished. This 
is qualitatively similar to the rapid evolution of a black 
hole interacting with a dipolar force- free configuration, as 
seen by [5T . In contrast with what is argued in [58 , mag- 
netic reconnection prevents the emergence of a relatively 
long-lived split-monopole configuration. In Sec. |VIII A| 
we discuss whether force-free or ideal MHD calculations 
more accurately describe the global structure of the cur- 
rent sheet, and the implications of the slower reconnec- 
tion that is seen in the MHD calculation of [57 . 

The evolution of the magnetic flux during the collapse 
is displayed in Fig. [Tt) again computed on a surface lo- 
cated at r = 1.5Rs. The total (signed) flux is small 
through most of the simulation, and only at late times 
does it becomes comparable to the absolute value of the 
magnetic flux. In the electrovacuum case, the final decay 
is governed by the main quasi-normal modes of the rotat- 
ing black hole. One observes faster decay in the force- free 



15 




-3-2-10123 -3-2-10123 

Ry Ry 



0.2- 




time ('ms'l 



FIG. 18: Rotating, unstable star. Time integral of the elec- 
tromagnetic luminosity, normalized to the peak EM energy 
of the magnetosphere, in both the force- free and electrovac- 
uum cases. The EM output depends weakly on rotation in 
the electrovacuum calculations, whereas in the force- free case 
the output is much larger 



FIG. 16: Rotating, unstable star (force- free) . Radial Poynting 
flux in red at t = (-0.3, -0.17, -0.01, 0.12) ms. The central 
colored zones mark the stellar fluid. The evolution of the 
poloidal magnetic field (blue lines) in the equatorial regions 
is qualitatively similar to that observed in the non-rotating 
case. 
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FIG. 17: Rotating unstable star (force- free) . The magnetic 
flux as a function of time, computed at r = 1.5i?s, normalized 
with respect to the initial value |$s(t = 0)|. The unsigned 
magnetic flux decreases after the black hole formation, similar 
to the non-rotating case. 
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FIG. 19: Rotating, unstable star (force-free). The elec- 
tromagnetic flux within annuli defined symmetrically in the 
northern and southern hemispheres by concentric cones hav- 
ing apertures in [ilS"", {i + 1)15''] {i = 0..5). The highest flux 
is within G [— 50^^, 50^^] . Here, as in the non-rotating case, 
the radiated energy decays faster than what would be ex- 
pected from a quasi-normal mode behavior, due to the effects 
of reconnection. 



run, just as as in the non-rotating case, possibly due to 
the greater self-inductance of the vacuum black hole. 
The radiated electromagnetic energy is displayed in 



Fig. 18, rescaled again with respect to the energy peak 
in the magnetosphere. 

The simulations indicate that around 20% of the en- 
ergy stored in the magnetosphere is radiated during the 
collapse in the force- free case, similarly to the electrovac- 
uum case. This is in clear contrast to the non-rotating 



case, where a much larger energy is radiated in the 
electrovacuum solution than in the force-free one. We 
note that although the energy of the magnetosphere is 
similar in the non-rotating and rotating cases (we find 
Cpeak = 1-5), the inclusion of rotation leads to a 20— fold 
enhancement in the EM energy radiated by a collapsing 
force-free magnetosphere: erad = 0.18. Hence Eq. (18) 



gives E^ad ^ 1.3 X 10^^ ^poie,i5 ^^S^ resulting in a strong 



average luminosity of L ^ 1.3 x lO^^^^^i^ erg s ^ dur- 
ing the collapse. 
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The distribution of the radiated energy is essentially 
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FIG. 20: Rotating, unstable star (electrovacuum). The elec- 
tromagnetic flux within annuli defined symmetrically in the 
northern and southern hemispheres by concentric cones hav- 
ing apertures in [il5°, (i + 1)15°] {i = 0..5). The radiated 
energy decays exponentially with a rate consistent with that 
expected from a quasi-normal mode behavior. 



quadrupolar, since the Newman- Penrose scalar ^2 has 
angular dependence mainly determined by an (/ = 2, m = 
0) mode. This radiation propagates outwards with a ve- 
locity V = 0.88c, which is very similar to what was ob- 
tained in the non-rotating case. As displayed in Fig. [19) 
most of the energy is radiated more efficiently near the 
angle = ±50^ with a peak intensity shortly after the 
formation of the black hole (For comparison purposes, 
Fig. [20] illustrates the behavior in the electrovacuum 
case). Subsequently this emitted radiation fades away 
rapidly as the magnetic flux is radiated away. 



VII. IMPLICATIONS FOR ASTROPHYSICAL 
TRANSIENTS: I. GAMMA-RAY BURSTS 

We now consider the implications of our simulation re- 
sults for high-energy transient phenomena. The very lu- 
minous gamma-ray bursts (GRBs) are generally believed 
to result from the formation of a stellar-mass black hole 
by sudden gravitational collapse, either in the core of a 
massive star [L,59,5 or after the merger of a neutron star 
with another compact object [2]|3]. A rapidly rotating 
magnetar is an interesting alternative [601462] . Our focus 
here is on ultraluminous EM outflows: the same general 
mechanism is strongly favored in other contexts such as 
pulsar synchrotron nebulae [63] and AGN jets [64 . 

Our calculations focus on the transition between a 
magnetar and a black hole, as a key part of the engine 
that drives a GRB. The implication is that a build-up 
of magnetic flux in the surface layers of a rapidly rotat- 
ing neutron star (formed, e.g., in a binary merger) can 
driven an electromagnetic outflow during the collapse of 



the star, independently of any Blandford-Znajek process 
operating afterward. As we explain here, this has some 
advantages over an EM wind operating before the col- 
lapse, in that this brief EM transient is likely to be sig- 
nificantly cleaner. We also compare the net EM output 
with that expected from a post-collapse jet, in the situa- 
tion where the surface layers of the neutron star, and not 
the surrounding disk, are the dominant zone for magnetic 
flux generated by a dynamo process. 

Our calculations have revealed strongly dissipative pro- 
cesses at work in the EM field surrounding the collaps- 
ing star. These include the reconnection of field lines 
near the magnetic equator, large-amplitude oscillations 
in the field, and (in vacuum calculations) the forma- 
tion of extended zones where :^ 5^, which in the 
presence of conducting matter imply relativistic motions 
of the entrained particles transverse to the magnetic 
field. The ejection of loops of magnetic field is observed 
in our non-rotating simulations. The dynamical evolu- 
tion of these loops has interesting implications for the 
gamma-ray flares from gravitationally stable magnetars 
(see Sec.lVmj). 



Although our calculations take into account the pres- 
ence of plasma implicitly by enforcing E • B = 0^ adding 
the associated force- free current and enforcing \E\ < \B\^ 
we expect that strong heating will occur in practice. 
Where the magnetic field reaches ~ 10^^ — 10^^ G, the 
implied temperature is above ^ 1 MeV. Such a high den- 
sity of electrons and positrons is created that macroscopic 
zones of non- vanishing E ■ B - diS are required by most 
pulsar discharge models (e.g. [65j) - cannot be main- 
tained. Effective heating can occur by other channels: 
for example, long- wavelength gradients in the magnetic 
field are converted efficiently to internal energy if the 
magnetic field becomes turbulent, so that a wide spec- 
trum of wave motions is formed that extends down to 
very small scales [66 . 

The relativistic outflow that is emitted by a rapidly 
collapsing star should therefore be quite hot. Large- 
scale magnetic fields and a relativistic photon-electron- 
positron plasma will both contribute substantially to the 
energy flux. Magnetar outbursts provide a fairly direct 
example of this phenomenon. 



A. EM Output Before, During, and After Collapse 

The star that collapses to form a black hole passes 
through three distinct phases: a pre-collapse phase dur- 
ing which it emits a magnetized wind if it rotates; the dy- 
namic collapse phase; and, if an orbiting disk is present 
- as it is following a neutron star merger - an accretion 
phase that is accompanied by a relativistic jet. Since our 
current simulations focus on the intermediate step, it is 
worth examining its relative contribution to the total EM 
output of the star. 

A large enough EM output A£^coiiapse is measured in 
our rotating force-free simulation to power some short 
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GRBs, especially if some account is made for beaming. 
Here the requirement is that the surface magnetic field is 
strong enough to hold off any accretion fiow. At an accre- 
tion rate M, this implies a polar magnetic field stronger 
thanill8 



pole 



MK(i^Ns) 



1/2 



7 X 10 



15 



M 



1/2 



Mo s-i 



G. (27) 



Here Vc is the circular speed, approximated as Keplerian. 
In this section we take a stellar mass Mns = 2.6 Mq 
and a radius i?NS = 15 km, as appropriate to a hot and 
rapidly rotating neutron star. The dipole field energy 
before the collapse is 



£^dipole,0 — Y2^Pole^NS 



1.5 X 10 



,49 



M 



M^s- 



erg. 



(28) 

We estimate AE'coiiapse ^ 0.3 ^dipoie,o, given that the 
energy released is about 0.2 times the peak magnetic en- 
ergy, which in turn is Cpeak ~ 1.5 times ^dipoie,o- From 
Eq. (I28|) we obtain 



collapse 



5 X 10"^^ 



M 



M^s- 



erg. 



(29) 



The precise numerical value depends non-linearly on the 
initial specific angular momentum J/Mns through the 
factor CradC'peak- faster initial spins imply stronger wind- 
ing of the magnetic field during the collapse. 

The magnetic field of an isolated neutron star acts as 
a couple between its reservoir of rotational energy and a 
dissipative outfiow. Even in the case of a (gravitationally 
stable) magnetar, the magnetic energy begins to domi- 
nate the rotational energy only at an advanced age, as 
the star spins down. Therefore the output of the pre- 
collapse phase could be very large compared with the 
release of EM energy during the collapse. Comparing 
the spindown energy radiated over a time At with the 
external magnetic energy (29) gives 



Lsd^t 

— =0.6 

^dipoie,o V ms y 



[Pns] ' f Rns V f At 



15 km 



NS 



(30) 



Here we have substituted the spindown power of an 
aligned, force-free rotator [67^ . 



^sd — 7^pole^NS<^ 



1 



(31) 



where (^ns = 27r/PNs- An isolated, rotating star would 
radiate energy equal to ^dipole, o in a few milliseconds.. 

In the context of binary neutron star mergers, one re- 
quires a magnetosphere to emerge from the very strong 



shear layer near the surface of the merger remnant (see 
Section IVH Bl for further discussion of how this could 
happen). If the neutron star is formed hot, a lengthy 
pre-collapse spindown phase would cause significant dif- 
ficulties with the application to short GRBs, because the 
wind generated during the pre-collapse phase is heav- 
ily loaded with nucleons and a particles that are driven 
outward by charged-current absorption of electron-type 
neutrinos near the neutrinosphere [54l [68l |69] . The con- 
nection between short GRBs and neutron star mergers 
would also be disfavored if this pre-collapse outfiow lasted 
longer than ~ 300 ms, given a characteristic short GRB 
lifetime of 0.03-0.3 seconds On the other hand, the 
survival of the merger remnant for ~ 100 — 300 ms would 
have the advantage of allowing stronger amplification of 
the magnetic field before the remnant collapses. The pos- 
sibility of longer-lived merger remnants for some configu- 
rations has been raised by recent simulations that employ 
a realistic, finite-temperature EOS [71 . 

Although two merging neutron stars are expected ini- 
tially to have magnet ospheres, the torus formed by the 
tidal disruption of the lighter star has a pressure vastly 
exceeding that of a typical pulsar dipole. The torus 
would, therefore, suppress a magnetosphere around the 
newly formed merger remnant. The neutrino-driven wind 
that emerges from the polar regions of the remnant will 
comb out the magnetic field, but this field need not ini- 
tially be coherent across the star or dynamically impor- 
tant. By the same token, if a torus were entirely absent, 
then the magnetic field threading the star would dissipate 
rapidly after the black hole forms, and the post-collapse 
phase would contribute negligibly to the output of the 
star. 

It is useful to express the EM power in terms of the 
"open" magnetic fiux that connects the surface of the 
star to the outfiowing wind. In the case of an isolated 
star, this is the fiux extending beyond the light cylinder. 



*^>open(^Ns) - TT^r (^^LC )^LC = *^NS 



(32) 



where <I>ns = ^poie • '^^ns dipolar magnetic fiux 



threading the star. One can then re- write Eq. (31) as 



1 



■^sd 



47r2c 



(^open^Ns) 



(33) 



After the star forms a distinct magnetosphere, the ac- 
cretion torus can continue to infiuence the wind power by 
modifying the open magnetic fiux. Let us suppose that 
the magnetic pressure dominates the torus ram pressure 
out to an equatorial distance Ra > Rns- Approximating 
the magnetosphere by a dipole, a fraction 



open 



NS 



Ra 



(34) 



of the stellar fiux is trapped by the torus. In this sit- 
uation, it is still possible for <^>open to exceed Eq. (32), 
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because the torus can extend inside the stehar hght cyhn- 
der, 



Rlc 

^open(no torus) Ra 



open 



3.3 



Ra 
15 km 



Pns 
ms 



(35) 

During our simulations of the collapse of an isolated 
star, we observe that the magnetic field lines are strongly 
twisted, so that most of the closed magnetic fiux opens 
out. Then 



(36) 



where Qh is the angular velocity of the horizon. 

A torus also plays an important role after the black 
hole forms by trapping a certain fraction of the mag- 
netic fiux that threads the star. If the magnetic field is 
strong enough to hold off the torus from the star before 
the collapse, then it will have a similar effect after the 
collapse. We therefore calculate the Blandford-Znajek 
power emerging from the horizon by assuming a uniform 
fiux density out to some radius i?^. From Eq. (8.65) of 
[72] . one gets. 
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(37) 



where (B) is the fiux density threading the region interior 
to the torus, and = '7t{B)R%. The difference in the 
normalizations of Eqs. (|33| and (37) largely refiects the 



fact that the torque on the black hole is maximized when 
the magnetic field has an angular velocity = VLh/"^' 
We can now relate Eq. ( 37 ) to the EM power generated 



before the collapse. The spin angular momentum is ap- 
proximately conserved during the collapse, J o:^ /ns^ns, 
and we also set Mbh = ^ns- Then 



J 



£l 



Pns 
ms 



Ml 



BH 



15 km/ \3Mq 

(38) 

where sj = /ns/^ns^ns ^ ^•3- The angular fre- 
quency of the black hole is ^hRh/c = J/MbylRsc = 
Jc/2GM^YL^ is related to I^ns pre-collapse by 
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NS 
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Ryl 
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(39) 



where Rsch = 2GMbh/c^ is the Schwarzschild radius. 
One expects (B) ~ 5poie after the collapse if most of 
the magnetic fiux threading the star is open before the 
collapse, and the pre-collapse magnetosphere is limited 
in size. The proportion of the trapped fiux threading the 
hole is 



H 



open 



H 



V^NS/ 



(40) 



We can now show that the EM power is suppressed im- 
mediately following the collapse. Substituting equations 
(|39|) and (pOl) into (IstI) gives 



0.01, 



(41) 



with a coefficient {Ra/ Rns)^ {Rns/ RaY{£i / 0.3)^ . Ad- 
ditional power will flow along the magnetic fleld lines that 
thread the ergosphere, but this portion of the black hole 
magnetosphere will mix with the accretion flow and may 
be less strongly magnetized. 

The interesting conclusion here is that the rapid twist- 
ing up of the magnetic fleld during the collapse can gen- 
erate a larger EM output than a fairly extended jet emis- 
sion after the collapse. It is worth summarizing the three 
main sources of this result: i) before the horizon forms, 
the EM power is proportional to rather than ^ ^^/4; 
ii) the magnetic flux remains pinned in the star for a few 
rotation periods during the collapse (before the onset of 
the black hole), and then springs out to flll a larger vol- 
ume; and iii) the relation between rotation frequency and 
angular momentum is enhanced by a factor ^ ej'^ ^3 
prior to the collapse; in other words, relativistic grav- 
ity has the effect of softening the growth of the rotation 
frequency as the star collapses. 

A reduction in the trapped flux (due to outward diffu- 
sion of the magnetic fleld into the torus) would, in this 
situation, initially increase the Blandford-Znajek power 
flowing from the horizon. The pressure of the trapped 
fleld approximately balances the ram pressure of the ac- 
cretion flow at some point outside the horizon. Therefore 
a reduction in the trapped flux allows the flow to reach 
closer to the black hole, and attain higher pressures, be- 
fore being interrupted. Although we are considering the 
flux originating in a dynamo process before the collapse 
(Sec. VII B), continued flux generation in the torus by 
the magnetorotational instability [73^ could play a role 
in modulating the jet power. 



1. Numerical Comparison 

Let us consider the EM energy that would be radi- 
ated following a binary NS merger, if the remnant sur- 
vives long enough to form a magnetosphere (Sec. VII B). 
Taking A£:coiiapse ^ 0.3E^ 

ipoie,0 5 and nor mali zing to 
an accretion rate O.IMq s~^ through Eq. (28), gives 
A£^coiiapse ~ 5 X 10^^ erg. We expect that the gain fac- 
tor A£^coiiapse/^dipoie,o depends non-linearly on the pre- 
collapse rotation rate, since the winding of the magne- 
tospheric fleld results from a competition between differ- 
ential rotation and torsional wave motion. In addition, 
A£^coiiapse depends indirectly on the torus mass and pres- 
sure through the strength of the magnetic fleld that is 
required to hold off the torus material from the neutron 
star surface. The larger accretion rate in a collapsar en- 
vironment implies a larger transient energy. 
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It is possible to make a direct comparison with the 
Blandford-Znajek jet that fohows the cohapse, if the 
magnetic flux that threads the black hole is left behind 
by the collapsing magnetar. Combining equations (30) 



and (41), and taking pre-collapse rotation period and 



radius 0.8 msec and 15 km, one finds that an energy 
~ 0.3£^dipoie,o would be radiated by a BZ jet over ~ 20 
ms. Equivalent ly, a Blandford-Znajek jet from a ~ 3 Mq 
black hole with Jc/GM'^^ ~ 0.7 would generate a power 
(|37l) Lbz ~ 1 X 10^^ erg/s. A recent binary merger 



simulation with GRMHD and dynamical gravity [74 pre- 
sented a polar magnetic field ~ 7 x 10^^ G developing 
from a much weaker seed field in the torus. The equa- 
tion of state used gave a fast collapse to a BH (within 
^ 10 ms), and therefore precluded the surface shear dy- 
namo that we have conjectured. (Note that this polar 
field could be affected by numerical resistivity, e.g. [75 ). 



2. Baryon Poisoning 



Comparing the output A£^coiiapse with Eq. (29), one 
sees that the pre-collapse star would release a compara- 
ble energy within ~ 100 ms after forming a magneto- 
sphere. However, it is well known (e.g. [54]) that such 
an outfiow would bear a much higher density of nucleons 
than a Blandford-Znajek jet from a black hole, due to the 
absorption of electron- type neutrinos and ant i- neutrinos. 
We expect that this nucleon loading would be strongly 
suppressed in the dynamical collapse phase, due to i) the 
relatively short duration of the emission; ii) the strongly 
wound field geometry {B^j^/Bp > 5); and iii) redshift- 
ing effects. Combining these effects suggests a significant 
suppression < 0.01 x 0.1 ~ 10~^ in the nucleon loading, 
so that the mass ejected would essentially be that present 
in the magnetosphere before the collapse. 



B. Emergence of a Magnetosphere via Dynamo 
Action in a Surface Shear Layer 

The immediate aftermath of a binary neutron star 
merger is distinguished from disk accretion onto a black 
hole, in that the velocity shear is strongest where the 
orbiting material makes a transition from centrifugal to 
hydrostatic support. Another feature which distinguishes 
the merger remnant from ordinary accreting neutron 
stars (X-ray pulsars) is that it does not initially have 
an ordered magnetosphere. 

Long after the first stage of the merger is complete, the 
velocity shear provides a tremendous source of free energy 
for amplifying a magnetic field. The power dissipated in 
the material settling onto the neutron star surface is 

1 



^shear ^ [V^ (R^s) - ^NS^Ns] 



3 X 10^^ 



M 



(42) 

where /shear = 1 - [^NS^Ns/Vc(^Ns)]^. Notc that the 
surface shear becomes more radially concentrated with 
time: as differential rotation is erased in the interior of 
the merger remnant, the surface shear is maintained by 
continuing accretion. The inner part of the shear layer 
develops positive dVt/dr and the magnetorotational in- 
stability is extinguished. 

A strong feedback mechanism is present which allows 
rapid magnetic field growth, but causes this growth to 
saturate once the star is able to form a magnetosphere 
that holds off the accretion fiow. When the magneto- 
sphere is present, the accreting material follows the mag- 
netic field and reaches the star at the same angular ve- 
locity. Shearing of the magnetic field in the outer layers 
of the star is therefore turned off. This effect is clearly 
demonstrated in the 3D accretion simulations of [82 . 

Only a tiny fraction of the accretion energy Lshear^^ 
must be converted to a poloidal magnetic field to hold 
off the accretion fiow: substituting (27) into the dipole 
energy (28) gives 
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Zl /shear erg S \ 



over a duration At. 

Low-mass X-ray binaries provide a nice example of sys- 
tems where this feedback appears to operate. The mil- 
lisecond radio pulsars that are descended from them have 
magnetic fields that are just strong enough {B ~ 10^ — 10^ 
G) to hold off an Eddington-level accretion fiow onto a 
neutron star - but not much stronger. The absence of 
persistent pulsations in the majority of LMXBs then re- 
quires that the magnetic field be aligned with the angular 
momentum of the accretion fiow. 

Now let us consider how the magnetic field evolves in 
the surface shear layer. The field present initially in the 
merging stars is rapidly amplified by a Kelvin-Helmholtz 
instability [28, or the magnetorotational instabil- 
ity [73]. Rapid growth of the magnetic field on large 
scales is sensitive to the speed of magnetic reconnection 
in the fiuid, and requires three-dimensional motions. The 
two-dimensional wrapping of a magnetic field by Kelvin- 
Helmholtz vortices does not generate net fiux; and the 
initial growth length of the MRI is very small compared 
with the scale height of the torus. MRI growth is fastest 
on a scale k'^ - B / {AirpY^'^Q. - 10-^^(5/10^2 G) r in a 
torus of mass - 0.01 Mq [73]. 

The hydrostatic structure of the inner shear layer 
makes it easier for the magnetic field to be pinned and 
retained than it would be in the surrounding torus. The 
magnetic field threading the shear layer is wound up, and 
since dVt/dr > 0, the mean toroidal flux density grows at 
least in a linear manner. If a hot, rotating, and massive 
neutron star (Mns ^ 2.6 — 3 Mq) can survive collapse for 
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more than ^ 100 ms, as is suggested by recent simula- 
tions of [TT, then this toroidal field becomes quite strong. 
Even if the seed poloidal field is as weak as 5p,o ~ 10^^ 
G - within the range of pulsar fields - then the toroidal 
field reaches 



[nc{RNs)-^Ns]AtBp^o 

At \ / Qc{Rns) — ^NS 



1 X 10^^ 



100 ms 
V 10^3 G 



10^ s- 



G. 



(44) 



Here CIc{Rns) is the angular frequency at the surface of 
the torus. 

Exponential magnetic field growth becomes possible 
when the wound-up field can rise buoyantly through 
the shear layer. The large dissipated power ( [42| ) will 
generate a strongly positive entropy gradient in the in- 
ner part of the shear layer. Where the material is 
convectively stable, only magnetic fields stronger than 

^ (GMNsMshear/i^^s)'/' ^ 1 X (^shear/O-l M©) ^2 

G can directly overcome the pressure of the over- 
lying material. But even in this case, the in- 
tense fiux of electron-type neutrinos allows a mag- 
netic field stronger than ~ 10^^ G to rise buoy- 
antly on the Alfven timescale, by erasing gradients in 
entropy and electron fraction that impede buoyancy 
[83j . The buoyancy time is tA ^ Ipl^npf^'^ / - 
1 {B(f)/l{)^^ G)~-^(p/10-^^ g cm~^)^/^ ms across a pres- 
sure scale height ~ i?Ns/4 ~ 3 — 4 km. Thus an 
exponential feedback loop appears quite likely in this sit- 
uation. 

One observes that the large-scale dipole magnetic field 
that emerges may be sensitive to the seed magnetic field, 
in the sense that a minimal seed field is required for 
the linearly wound field to reach the buoyancy threshold. 
Nonetheless, it is also clear that magnetar-strength mag- 
netic fields do not require magnetar-strength seed fields 
in the presence of persistent surface shear. 

The strength of the magnetosphere that eventually 
emerges results from a competition between the finite 
rate of amplification and the diminishing shear stress at 
the surface of the star. At an age of ~ 10(100) ms, a 
poloidal magnetic field stronger than ~ 10^^ G (10^^ 
G) will begin to apply a strong, negative feedback on 
the surface shear. Once the torus mass drops to the 
point that accretion is mediated mainly by magnetic 
stresses, the accretion rate is Mq ^ ^T,o/^diff,o ^ 
10(a/0.1)(MT,o/0.1 Mq)Mq s-^ Here ' Mt, o ' and 
^diff,o ^ 10(a/0.1)~^ ms are the initial torus mass and 
diffusion time in this viscous phase, and a is the viscosity 
coefficient. The accretion rate drops as the torus mate- 
rial spreads outward, as M ~ ^o(V^diff,o)~^^^, as long 
as the torus remains geometrically thick and conserves 
its angular momentum (e.g. f84j)- 



C. Galactic Magnetars and Delayed Collapse 

We now consider the observational imprint of a magne- 
tar if it forms by the accretion of a thin layer of strongly 
sheared material, but is (initially) gravitationally stable. 
We focus on the rate of energy loss and the spin history 
of the star. 

There is a substantial reduction in the energy release 
by spindown, compared with a star that initially ro- 
tated as a solid body, due to the internal rearrange- 
ment of angular momentum. Approximate solid-body 
rotation is attained on the poloidal Alfven timescale 
i?Ns(47rp)^/^/5poie, which is generally much shorter than 
the magnetic-dipole spindown time 
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Most of the shear energy is then dissipated internally. 

A delayed collapse is possible if the magnetar only 
slightly exceeds the maximum mass for a non-rotating, 
zero-temperature neutron star. The loss of rotational 
support by a magnetic wind would trigger a collapse on 
the rotational braking time. 

Quite generally, the accretion of a thin layer of strongly 
sheared material provides an attractive mechanism for 
creating magnetars, and so the results of this section 
should have a broader application to the Galactic mag- 
netar population, even if these stars are well below the 
maximum mass. Several possible channels are available: 
the rotation of the accreting material could be gener- 
ated by an instability of the supernova shock [85 -88 ; or 
the neutron star could be exposed to material with very 
strong vorticity in a merger event, e.g. when it merges 
with a companion CO white dwarf, or with the core of 
an evolving Be star. 

The spindown power that is released by the magnetar 
is then reduced significantly compared with a star which 
rotates uniformly with the same angular velocity as the 
surface material. The final spin period resulting from the 
accretion of a layer of mass Mghear and rotation period 
^shear outo a ucutrou star of total mass Mns is 
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where jshear = (Mshear/O.lM0)(Pshear/ms) ^ As the 

star spins down, it deposits a rotational energy 

(47) 

in the surrounding shock wave. For a magnetar of po- 
lar field Ppoie, the corresponding spindown luminosity is 
fairly modest. 
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and the star enters the extended spindown phase after a 
time 
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In this situation, only a very strong internal magnetic 
field (^toroidal ^ 10^^ G) would induce sufficient triaxi- 
ality in the star that the gravity wave torque competed 
with the external electromagnetic torque |89] . 



VIII. PLASMOID DYNAMICS AND 
RADIATION: IMPLICATIONS FOR MAGNETAR 
FLARES 

Escaping loops of magnetic field are observed in our 
non-rotating simulations. These loops are formed by 
magnetic reconnection, which means their size and rate 
of formation are sensitive to the treatment of resistiv- 
ity, and to the disregard of stresses associated with bulk 
plasma fiow by the adoption of the force- free equations. 
An outgoing plasmoid (containing plasma and a closed 
magnetic field) has a well-defined center of mass: we find 
that its bulk motion is measurably smaller than the speed 
of light, Vbuik ^ 0.9c and Fbuik ~ 2, at a distance r ~ 10 
times the pre-collapse neutron star radius R^s- This 
is significantly less relativistic than the bulk motion ex- 
pected for a gas of freely expanding particles released by 
the collapsed star (Fbuik ^ ^/Rns ^ 10). 

After addressing each of these issues, we make contact 
with the giant gamma-ray fiares of the Galactic magne- 
tars. These appear to involve the ejection of an energetic 
plasmoid (~ 10^^ — 10^^ erg), but due to the shearing of 
the magnetic field lines rather than the collapse of the 
star 190J. 



A. Dependence of Reconnection Rate on 
Resistivity Model 

When magnetic field lines of an opposing sense are 
stretched out and forced into contact, the rate at which 
they reconnect is sensitive to the microscopic model of re- 
sistivity. In an ohmic plasma with a uniform resistivity, 
a long current sheet forms and reconnection is very slow; 
fast reconnection with an x-point geometry depends on a 
local maximum in the resistivity [91 . In some contexts, 
such as the Solar corona, the microscopic explanation for 
this behavior may be provided by the Hall terms in the 
conductivity [92 . These are relatively less important if 
the plasma is loaded with e+/e~ pairs, as would be ex- 
pected in the magnetosphere of a collapsing magnetar. In 



a fiuid, small-scale hydromagnetic turbulence appears to 
greatly accelerate the reconnection rate [93 ; but whether 
such a process can operate in a low-/3 plasma of astro- 
physical dimensions is not yet determined. 

In the present work we find evidence for relatively fast 
reconnection of magnetic field lines, V^ec ^ 0.1 c, as 



calculated in the force- free approximation (Sec. VIA). 
This results from a change in topology of the field lines 
(the formation of an x-point). It is not due to strongly 
enhanced dissipation in an extended equatorial current 
sheet: E.J dissipation is measured to be small in both the 
rotating and non-rotating cases. Such a concern arises in 
force- free calculations of pulsar magnetospheres, where 
the magnetic field is less dynamic and is anchored in 
the star. There the absence of fiuid pressure support 
in the current sheet leads to a rapid collapse of the mag- 
netic field toward the sheet, unless explicitly compen- 
sated [9]. In the present case, after the black hole 
forms the infiow of magnetic fiux toward the equator can 
continue in a more dynamic manner into the star, or out 
to the computational boundary. 

A recent treatment of reconnection in the magneto- 
sphere of a stationary black hole [57] illustrates a slower 
field decay when employing an ideal MHD treatment, as 
compared to a force-free approach. (Resistivity in the 
MHD case arises through the numerical approximation.) 
Such a slower decay is also observed at late stages after 
the formation of the black hole in our simulations when 



comparing force-free and ideal MHD fields |1 19] : though 
by this time the field strengths are orders of magnitude 
below their peak values. A comparison of force-free and 
MHD reconnection calculations which explores more gen- 
eral resistivity models, and their infiuence on the recon- 
nection geometry, remains to be developed. 



Magnetic Reconnection Delayed 
by Plasma Outflow 



Reconnection is usually studied in a context where the 
plasma fiow speed along the magnetic field is a small frac- 
tion of the Alfven speed: for example, in a steady MHD 
wind, the conservation of angular momentum implies a 
slow drift of particles along the spiral magnetic field out- 
side the Alfven critical point. However, in some contexts, 
such as magnetar fiares, the fiow speed can approach the 
speed of light. Similarly, in our collapse simulations we 
see large-amplitude motions on the magnetic field loops 
threading the neutron star, which suggest strong plasma 
heating. The backreaction of the outfiowing plasma on 
the reconnection of field lines is not taken into account. 

Even when the magnetic field has a tendency to re- 
connect through an x-point, reconnection will be delayed 
until the kinetic pressure ^ of the outfiowing pair- 
photon plasma (with thermal energy density ~ U) drops 
below the Foynting fiux that would fiow toward the cur- 



22 



rent sheet in the absence of plasma flow. One requires 



C//?2 < 0.1- 



c Stt' 



(50) 



where the coeflicient on the right- hand- side is appropri- 
ate for fast x-point reconnection. 

The magnetic fleld hnes are stretched outward by the 
plasma flow beyond an Alfven radius |66] 
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given i?NS ^ 10 km. Note that, at the peak of the out- 
flow, the pressure of the pair-photon fluid is comparable 
to the pressure of the stretched magnetic fleld lines at 
r = Ra^ and outside this radius grows as ~ {v/RaY 
with respect to the split-mo nopole fleld pressure. For ex- 
ample, in a magnet ar giant flare, the ~ 0.1 s width of 
the main gamma-ray pulse is comparable to the time for 
magnetic and elastic stresses to rearrange the stellar in- 
terior; but it is ~ 300 times longer than the flow time 
out to the Alfven radius (51), and therefore much longer 
than the timescale for x-point reconnection at a speed 
- 0.1c. 



C. Radio Afterglow from Strongly Magnetized 
Outflows 

A magnetically-dominated plasma that is ejected from 
a collapsing magnet ar (or a nearby Soft Gamma Re- 
peater) can be a much stronger source of synchrotron 
emission than a shocked plasma of comparable energy 
density. Existing calculations of radio afterglows of short 
GRBs ([94l and references therein), as well as calcula- 
tions of the radio afterglow of SGR giant flares [95 , focus 
on synchrotron emission by a population of non-thermal 
electrons that are accelerated at a shock that leads the 
outflow. For the bulk of GRBs, one infers efliciencies of 
conversion £e^£B ^ 0.1 of bulk kinetic energy to non- 
thermal electrons and to magnetic flelds downstream of 
the forward shock. These moderate values of £e and £b 
are consistent with the broadband tails of radio-X-ray 
emission that follow the brief, bright gamma-ray phase. 

The synchrotron emission from a magnetically- 
dominated plasmoid will be proportionately much 
brighter, by up to a factor ~ 100, for two reasons. First, 
and most obviously, sb now approaches unity. Second, if 
the magnetic energy density dominates the thermal en- 
ergy density, then a cascade process (involving the cre- 
ation of high-wavenumber Alfven modes) preferentially 
heats the electrons [96 . If the plasma is very relativistic, 
damping is mainly due to charge-starvation of the waves, 
at wavenumbers where the amplitude of the fluctuating 



current density begins to exceed erieC [66[ |77]. (A differ- 
ent damping mechanism operates at higher ion densities: 
the waves are Landau-damped on the parallel motion of 
the electrons.) In practice, the relative amplitudes of the 
bulk synchrotron emission and the shock emission will 
depend on the degree of disorder in the plasmoid mag- 
netic fleld. But for magnetar flares, indirect evidence 
that bulk synchrotron emission dominates comes from a 
rapid initial drop in radio flux that is consistent with the 
sudden compression of the plasmoid, followed by rapid 
adiabatic cooling [95^ . 



D. Applications to Magnetar Outbursts 

Magnetar outbursts involve more limited releases of 
energy that leave the original star intact. They are trig- 
gered when the footpoints of a ~ 10^^ G magnetic fleld 
are strongly sheared by an elastic instability of a neu- 
tron star crust, thereby generating a hot plasma and 
an intense burst of gamma-rays [90^ 97 . The flrst, ex- 
tremely bright, stage of a giant magnetar flare lasts only 
~ 100 ms and bears a considerable resemblance to a 
'classical' gamma-ray burst, demonstrating high temper- 
atures {kT > 200 keV) and a signiflcant non- thermal 
component to the spectrum [98 . (This phase saturates 
almost all X-ray detectors, but was well-resolved by the 
Geotail experiment |99J in the 27 August 1998 and 27 
December 2004 flares.) The duration and luminosity are 
consistent with the internal rearrangement of the mag- 
netic fleld in a neutron star, with a strength ~ 4 — 5 x 10^^ 
G based on considerations of magnetic fleld transport and 
global flare energetics [100 , several times stronger than 
the standard dipole expression for the spindown-derived 
magnetic fleld. 

The giant flares appear to involve the ejection of a plas- 
moid. The combination of fast variability and extreme 
luminosity (up to ^ 10^^ ergs s~^) implies, through the 
usual arguments of gamma-ray opacity |101j , that the 
emitting plasma has expanded to a much larger volume 
than that of the neutron star. This expanding plasma is 
an excellent electrical conductor and must carry some of 
the stellar magnetic fleld with it. 

This expected property of magnetar giant flares helps 
to explain two apparently contradictory phenomena. A 
straightforward argument based on the theory of thermal 
flreballs (e.g. |102j ) shows that the expanding plasma is 
moving with a high Lorentz factor at the radius (^10^ 
cm) where it becomes transparent to the gamma-rays 
j97] . The observation of variability on a timescale 5tvar ^ 
4—20 ms in the gamma-ray flux implies strong constraints 
on the baryon rest energy flux. The advected electrons 
and ions must become transparent to the gamma-rays 
close enough to the magnetar that the differential light- 
travel time r/2r^c across the outflowing plasma is shorter 
than ^tvar- Given a total outflow luminosity L, this im- 
plies that the ion rest energy contributes no more than a 
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fraction 

L \ Lgt J 

of the energy flux. 

On the other hand, radio monitoring detected transient 
emission in the weeks fohowing both flares |1Q3[ llQ4j . 
The emission following the 2004 flare was especially 
bright, as befitting the much greater energy of the burst, 
and could be tracked on the sky for more than a year. 
The expansion of the radio source implies a transverse 
velocity v ~ 0.7(1^/15 kpc)c [104 . After this, clear evi- 
dence is seen for a break in the radio light curve consistent 
with a transition from uniform expansion to a Sedov- 
like phase. We infer that the measured proper motion 
is probably the free expansion velocity, as corrected for 
relativistic aberration. 

In principle it is possible for the measured transverse 
speed to be less than the speed of light, if the intrinsic 
motion V is nearly luminal but the motion is directed 
away from the observer (at some angle 6 > it/2 with 
respect to the line of sight to the magnetar): /3^,obs = 
/3sin^/(l — j3 cos 0). The observed relative brightness of 
both the radio and gamma-ray emission argues against 
this: one measures ^ A x 10'^^(i:>/15 kpc)^ in 2004 
vs. r^bx 10'^^(D/15 kpc)2 in 1998, and a radio flux 
Fj^ ^ 50 mJy vs. 0.3 /iJy at 8.5 GHz 1 week after the 
flare [Ml Il03l Hog. (A factor - 2 underestimate of the 
~ 15 kpc distance to SGR 1806-20 appears unlikely given 
that the source position is in the Galactic plane.) 

The energy reservoir that powers the early, subluminal 
stage of the radio afterglow must be composed of some- 
thing other than electrons and positrons emitted by the 
star, which would mainly have annihilated during the 
very brief fireball phase. A reconnected magnetic field 
is the most plausible delayed carrier of energy, especially 
given the strong limitations on the baryon flux during 
the < 100 ms gamma-ray pulse. A closed loop of mag- 
netic field carries a finite inertia, and so its center-of-mass 
frame will not accelerate with distance from the source 
as does that of a collimated particle beam. 

The outflowing pair-photon fluid must overcome the 
tension of the magnetic field lines that are anchored in 
the star, and so comparable energy can be put into the 
stretched field, which is pulled into a split-monopole con- 
figuration. The energy of the stretched field is concen- 
trated close to the star, B'^r^ ~ r~^, although not as 
strongly concentrated as in a static dipole {B'^r^ ~ r~^). 
In a static plasma, the time for reconnection at the Alfven 
radius Ra is tree ~ 3(7^^/100 km)(VA/0.1 c)~^ ms. 
Comparing this expression with the measured e-folding 
time of ~ 30 ms for the tails of the giant flare pulses 
[99] suggests that reconnection is gated by the decrease 
in pair-photon pressure. 

We conclude that, in a magnetar giant flare, most of 



the released magnetic energy tends to follow the pair- 
photon pulse, with the delay between the two compo- 
nents being due to continued shearing of the external 
magnetic as the interior of the magnetar adjusts on the 
internal Alfven time of ^ 0.1 s. 



IX. CONCLUDING REMARKS 

Understanding the global behavior of strongly gravi- 
tating, dynamical systems, containing dense matter cou- 
pled to ultra-strong magnetic fields, is of key importance 
for a thorough understanding of possible signals from 
them. In this work, we have presented a new approach 
to this end by combining the ideal MHD and force-free 
approximations in a suitable manner within general rel- 
ativity. 

We expect that this hybrid scheme represents real 
progress towards greater realism. The stellar interior 
utilizes ideal MHD to faithfully model the neutron star 
without the disadvantages incumbent in the low density 
exterior. Likewise, the exterior solution uses the force- 
free approach and so captures the dynamics of the ten- 
uous plasma. The entire domain is described by a fully 
nonlinear and fully dynamic general relativity solution 
necessary for strong-field gravity. Indeed, a key aspect 
of this approach is its generality. It does not require a 
prescribed stationary stellar boundary, and can be ap- 
plied, for example, to dynamical systems such as collaps- 
ing stars and non- vacuum compact binaries. 

We have exploited this approach to study stellar col- 
lapse in both rotating and non-rotating cases. Our stud- 
ies reveal a rich phenomenology in the magnetosphere as 
the collapse proceeds. In particular, magnetic reconnec- 
tion plays an important role by inducing strong electro- 
magnetic emission as well as the infall of electromagnetic 
energy into the black hole, which in a short time loses 
all its 'hair'. When the star starts off rapidly rotating, 
the energy of the magnetospheric plasma grows signif- 
icantly during the collapse. As the star rotates faster, 
its magnetic field lines do not have time to re- adjust to 
the increased rotation, and are strongly wound up out 
to the initial light cylinder. It is worth emphasizing that 
this conversion of dynamical energy into electromagnetic 
energy does not depend on resistive effects. 

Two issues of principle have arisen in performing these 
calculations. First, we have shown conclusively that the 
force- free approximation to the evolution of a dilute, rel- 
ativistic plasma inevitably leads to singularities. These 
singularities are avoided in our calculations by continu- 
ously pruning the electric field. In the simulations that 
we have run, this procedure appears to cause limited en- 
ergy dissipation, but its necessity should be kept in mind. 

The second issue of principle regards the dependence 
of the reconnection geometry on the resistivity model. 
Fluid pressure is responsible for slowing down the rate of 
reconnection unless an x-point geometry is able to form. 
By neglecting fluid pressure, the force-free approxima- 
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tion clearly facilitates the formation of x-points. In spite 
of this, we observe very limited numerical dissipation in 
current sheets. Although the plasma that is represented 
by our model magnetospheres is, in reality, strongly colli- 
sional, it should be kept in mind that fast x-point recon- 
nection still occurs in collisional plasmas in the presence 
of hydromagnetic turbulence - as has been demonstrated 
so far in weakly magnetized plasmas [93 . Some further 
exploration of the reconnection geometry is possible in 
resistive MHD calculations of stellar collapse by varying 
the spatial dependence of the resistivity. 

We have also discussed how the electromagnetic out- 
bursts from collapsing magnetars may have interesting 
observational effects. In particular, our results are rele- 
vant to binary neutron star merger scenarios, in which 
a hypermassive neutron star forms and collapses to a 
black hole. Collapse to a black hole can happen either 
promptly or after many dynamical times, depending on 
the masses involved and the equation of state describ- 
ing the stars (see, for instance, [19]). Although Kelvin- 
Helmholtz and magnetorotational instabilities will create 
ultrastrong magnetic fields io^^~^^ G), we have sug- 
gested that the global field (in particular, the amount of 
magnetic fiux threading the merger remnant and even- 
tual black hole) could depend strongly on the lifetime of 
the remnant. The formation of a distinct magnetosphere 
would require tapping only ~ 10~^ — 10~^ of the en- 
ergy dissipated in shear layer at the remnant surface. As 
such a star collapsed, the magnetosphere would qualita- 
tively follow the behavior outlined here. The electromag- 
netic output of the collapse could compete with the later 
emission from a Blandford-Znajek jet emanating from the 
black hole horizon, and source a powerful electromagnetic 
counterpart to the gravity wave signal (e.g. |lQ6j ). In ad- 
dition, the magnetic field dynamics that is revealed in our 
simulations has a number of interesting implications for 
gamma-ray bursts and magnet ar fiares, as discussed in 
Sees. EUl and Eini 

Beyond the work analyzed here, our approach is readily 
applicable to other relevant systems and will be applied, 
in particular, to study binary neutron star systems ^lQ7j . 
and black hole- neutron star systems (see e.g. |lQ8j ). It is 
important to stress, however, that our approach is not 
free of ambiguities, in particular with respect to how 
and where the matching between the force-free and ideal 
MHD regions is implemented. To the extent possible, we 
have tested the robustness of our results versus known 
solutions, which give us confidence in this approach. 
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Appendix A: The Transition from Ideal to 
Force-Free MHD 

Details of the transition from the ideal MHD regime 
to the force-free regime merit particular attention, espe- 
cially in our rotating collapse solutions. Since we match 
two different formulations of Maxwell's equations coupled 
to conducting matter, it is important to understand how 
the choice of matching layer affects the exterior force- free 
solution. Of course, an unambiguous test can only be 
provided by a complete resistive MHD solution that can 
handle the strongly magnetized regions outside the star 
and, in particular, can follow the large changes in density 
and rotation that are encountered during the collapse. 
While work on this direction is in progress, the simpler 
approach presented here allows us to obtain a first so- 
lution to the relativistic magnetosphere in this strongly 
dynamic situation. 

The positioning of the matching zone is constrained by 
competing considerations. On the one hand, if it sits too 
close to the surface, then the interior MHD solution that 
sources the exterior force-free solution will be unrealistic, 
due to the density fioor that is applied in the MHD at- 
mosphere. One thus might want to place the transition 
layer well within the stellar surface. However, such a deep 
layer might underestimate the magnetic field strength at 
the base of the magnetosphere, and imply force- free be- 
havior of the magnetic field where that approximation is 
not justified 

In particular, we have found that placing the transi- 
tion zone at too low a density implies an unrealistically 
large toroidal magnetic field at the base of the force-free 
zone. Recall that the region exterior to the star, with 
its tenuous plasma, will have a relatively large magne- 
tization, and should be forced to co-rotate with the in- 
terior of the star. Unless prohibitively high resolution 
is employed, the atmosphere of the MHD solution has a 
large enough inertia that this condition can be violated. 
Instead, the magnetic field experiences a non-negligible 
(and non-physical) differential rotation. This effect is 
easily exacerbated in a collapsing solution. 

Within our current approach, we address these issues 
in two ways, by (i) adopting an initial configuration of 
the force-free fields and ideal MHD fields respecting such 
co-rotation up to the light cylinder and with a suitable 
differential rotation outside that radius; and (ii) choosing 
a transition layer at an appropriate distance inside the 
stellar surface. 

Regarding the initial data, since the rotation induces 
a toroidal magnetic field we adopt initial data that has a 
toroidal magnetic defined as 5^ = —QpRcyiBr and also 
E = —V X B throughout domain. We also choose the 
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lines so that = l^star near its surface, propagating 
along the magnetic field lines with a speed ^ c. 

The sharp negative gradient appearing in the last 
stages of the collapse refiects the strong radial gradient 
in l^star that appears near the rotation axis, as well as 
the onset of strong general relativistic effects. As the 
event horizon arises (prior to the formation of the appar- 
ent horizon at t = 0), it disconnects the interior from the 
exterior solution and causes l^i? to decrease, tending to 



FIG. 21: Rotating, unstable star (force- free) . Isosurfaces 
of the kernel function of Eq. ([9| at two different times t = 
—0.30ms, left, and t = —0.07ms, right, while the rotating 
star collapses. Mapped in color is the density of the star in 
cgs units in a eridional plane. Black lines mark two isosurfaces 
of the kernel function, corresponding to F = (0.01, 0.99). 



MHD atmosphere to co-rotate with the star for Rcy\ < 
2i^star, and not to rotate otherwise. 

Both the location of the stellar surface and the density 
in the transition layer are dynamical. It is important 
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that the choice of transition density 
spects the increasing density of the collapsing star. In 
order to maintain a consistent depth of the transition 
layer, its position is adjusted as the star collapses by 
scaling the matching density p^^^'^^ in proportion to the 



peak density within the star: p' 
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0)[pmax(^)/Pmax(^ = 0)]. The transition layer is dis- 
played at the beginning and near the end of the collapse 
in Fig. |2l] Notice that with this conservative approach 
we are underestimating possible rotational effects. 

The angular frequency ftp of magnetic field lines an- 
chored near the rotation axis of the star is illustrated in 
Fig. [22] From the initial condition (in which constancy 
of l^i? is enforced), a negative gradient in l^i? develops 
in the transition layer during the earliest stages of the 
collapse. This negative gradient has then disappeared by 
the time that the stellar angular velocity has increased 
by ~ 20%. above the initial value (see Fig. 14). 
From then on, the qualitative behavior obtained is con- 
sistent with the expected one: during the collapse, the 
star transfers angular momentum to the magnetic field 
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FIG. 22: Rotating, unstable star (force- free) . (Top panel) 
Angular velocity Qf oi the magnetic field line emanating from 
= 122°, = 0° as a function of cylindrical radius (normal- 
ized with respect to its initial value). (Bottom panel) Shape 
of such magnetic field line as time progresses. 



the value = l^i^/2 expected by the Blandford-Znajek 
solution of a spinning black hole. 
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